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Advanced  hUersonictaii-breathing  and  hybrid  engines  utilizing  subsonic/supersonic 
combustion/,  and  hlgtijmergy  lasers,  are  required  for  effective  weapon  systems  in 
future  aertospace/Aissions.  The  overall  objective  of  this  prograa  is  to  investigate 
the  turbulent  raixfftg  and  combustion  processes  required  for  these  advanced  systems. 
Results  oft  these  studies  will  assist  in  engineering  design  and  development  o£ 
practical  injection-mixing  systems  and  combustion  chambers  for  the  advanced propulsi — 
systems,  as  well  as  providing  critical  design  input  for  development  and  optimiza¬ 
tion  of  high  energy  chemical  lasers.  The  following  theoretical  ane  experimental 
: ud.es  .cr  homogeneous  end  heterogeneous  highspeed  turbulent  mixing  and  combustion  i~ 
subsonic  and  supersonic  streams  are  reported.  1)  A-'  analytical  study  of  turbulent 
reacting  flows  with  mass,  momentum,  and  energy  transfer,  2)  Experimental  and 

•S'***'  C^-particlc,  turbulent  mixing  and  reacting  flow  with  non-tairger.- 
P  f  ^ ^  J  -  on.  3)  Analysis  of  injection  of  non-homcgar.eous  and  ha  ter  a- 
g  ous  fuels-  into  unnorm  supersonic  and  subsonic  streams,  and  4)  A  review  of  exxstir 
experimental  data  and  analytical  methods,  concerning  turbulent  mixing  and  co=busci~- 

Objectives  of  this  research  are  to:  1)  Determine  effects  of  chemic-l''re’c*ioa  _ 

momentum  end  energy  transport  on  flow  field  development  of  coaxial  jets  and  ducted  ' 
°MelOP  ,iecessa^  analytical  tools  for  predicting  gas-particle  nixing  and 

Imsm ^and^r^i  3)  Inve®Ll8ate  s«ling  parameters,  pressvre  gradient  shock 

losses,  and  Initial  i*t  penetration  versus  downs  trea-i  turbulent  diffusion  for 
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ABSTRACT 


Advanced  hypersonic  air-breathing  and  hyt  rid  engines  utilizing  subsonic /supersonic  combustion, 
and  high  energy  lasers,  are  required  lor  effective  weapon  systems  in  future  aerospace  missions.  The  over¬ 
all  objective  of  this  program  is  to  investigate  the  turbulent  mixing  and  combustion  processes  required  lor 
these  advanced  systems.  Results  of  these  studies  will  assist  in  engineering  design  and  development  of 
practical  iniection-mixing  systems  and  combustion  chambers  for  the  advanced  propulsion  systems,  as  well 
as  providing  critical  design  input  for  development  and  optimization  of  high  energy  chemical  lasers. 

The  following  theoretical  and  experimental  studies  for  homogeneous  and  heterogeneous  high¬ 
speed  turbulent  mixing  and  combustion  in  subsonic  and  supersonic  streams  are  reported: 

( I  )  An  analytical  study  of  turbulent  reacting  flows  with  mass,  momentum,  and  energy  transfer. 

(2)  Experimental  and  analytical  study  of  gas-particle,  turbulent  mixing  and  reacting  flow  with 
non-tangential  particle  injection, 

f  3 1  Analysis  of  injection  of  non-homogeneous  and  heterogeneous  fuels  into  uniform  supersonic 
and  subsonic  streams,  and 

(  4)  A  review  of  existing  experimental  data  and  analytical  methods,  concerning  turbulent  mixing 
and  combustion. 

Objectives  of  this  research  are  to: 

{ I )  Determine  effects  of  chemical  reaction,  mass,  momentum,  and  energy  transport  on  flow 
field  development  of  coaxial  jets  and  ducted  flows, 

(2)  Develop  necessary  analytical  tools  for  predicting  gas-particle  mixing  and  combusting  flow 
systems. 

f3)  Investigate  scaling  parameters,  pressure  gradient  shock  losses,  and  initial  jet  penetration 
versus  downstream  turbulent  diffusion  for  the  transverse  mode  of  injection,  and 

(4)  Determine  the  adequacy  and  applicability  of  existing  two-stream  mixing  data  and  anah  tical 
techniques. 
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I.  SUMMARY 


The  overall  objective  of  this  program  is  to  provide  the  flnn  technological  basis  required  for  the. 
design  and  optimization  ot  advanced  propulsion  systems  and  high  energy  chemical  and  gas  dynamic  lasers, 
A  dual  approach  is  being  followed  to  achieve  this  objective:  I )  the  gaseous  turbulent  nixing  and  com¬ 
bustion  processes  are  being  characterized  for  flow  conditions  and  geometries  of  practical  interest,  and 
2)  the  injection,  mixing  (dispersion)  and  combustion  processes,  required  for  the  effective  utilization 
of  high  performance  powdered  fuels,  also  arc  being  investigated. 

A.  TURBULENT  MIXING  AND  REACTING  FLOW  CHARACTERIZATION 

In  diffusion  flames  the  mixing  or  transport  of  mass,  momentum,  and  energy,  and  chemical 
reaction  occur  simultaneously  during  the  combustion  processes.  These  flames  are  generally  turbulent 
because  of  the  significant  shear,  and  hence  turbulence,  generated  during  the  injection  and  combustion 
processes.  Such  flows  cannot  be  characterized  without  experimental  data,  which  are  required  to  specify 
the  Reynolds  transport  terms  for  each  specific  class  of  flows,  e.g.,  free  jets.  A  general  semiemptrical 
approach  is  presented  for  modeling  these  terms  in  practical  turbulent  reacting  and  nonreacting  flows. 

Results  are  presented  which  demonstrate  that  the  transfer  of  mass  (mass  mixing)  is  significantly 
more  important  than  the  transfer  of  either  energy  or  momentum  for  a  coaxial  hydrogen  jet  reacting  with 
an  external  high-speed  air  stream  These  results  demonstrate  that  the  key  to  successful  prediction  of 
practical  combustor  and  chemical  laser  performance, when  diffusive  burning  is  employed,  is  the  determin¬ 
ation  of  suitable  mass  transfer  models,  rather  than  momentum  transfer  models  e.g.,  for  the  eddy  viscosity. 
The  need  to  specify  a  variety  of  turbulent  mass  transfer  coefficients  for  various  chemical  species  for 
characterization  of  simultaneously  mixing  and  reacting  flows  was  suggested  by  a  fundamental  argument 
based  on  a  general  consideration  of  the  interrelation  of  turbulence  and  chemical  kinetics,  and  more 
directly  by  an  analysis  of  gas  sampling  data  reported  for  a  specific  hydrogen-air  reacting  flow. 

The  technique  for  direct  determination  of  turbulent  transport  coefficients  by  differentiation  of 
experimental  data,  designated  the  Inverse  Solution  Technique,  was  applied  to  a  number  of  nonreacting 
flows.  Results  demonstrated  that  valid  transport  coefficients  could  be  obtained  for  all  of  these  data, 
and  that  the  Composite  Trends  resulting  from  this  analysis  are  very  useful  in  the  development  of  practical 
turbulent  mixing  models.  This  technique  was  successfully  extended  to  analysis  of  a  turbulent  reacting 
hydrogen-aii  flow.  The  resulting  hydrogen  mass  transfer  coefficient  was  shown  to  be  very  adequate  for 
prediction  (numerically )  of  the  experimental  hydrogen  concentration  profiles  used  in  its  determination. 

In  addition,  these  results  were  utilized  to  modify  empirically  an  existing  eddy  viscosity  model  developed 
for  nonreacting  flows,  so  that  adequate  agreement  between  predicted  and  experimental  profiles  could  be 
attained  for  the  reacting  flow.  Finally,  .esults  obtained  using  the  Inverse  Solution  Technique  were  shown 
to  be  useful  for  the  se*cction  of  the  most  satisfactory  mixing  model  applicable  to  a  particular  flow  condi¬ 
tion  and  geometry. 

B.  MODELING  THE  (  ORE  REGION 

Two  eddy  viscosity  models  applicable  to  the  core  region  of  a  turbulent  coaxial  jet  were  developed 
by  modifying  existing  transition  region  models.  This  region  is  of  critical  interest,  because  in  diffusion 
flames  it  is  generally  in  this  region  that  the  ignition  and  most  of  the  combustion  occurs.  Each  of  the 
models  was  demonstrated  to  apply  to  select  data  when  either  simple  step-type  {slug)  initial  profiles 
were  used,  or  when  experimental  profiles  in  the  core  were  used  to  initiate  the  numerical  integration. 
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These  models  were  successfully  applied  over  a  ranee  of  flow  conditions  ami  geometries;  however,  further 
effort  is  required  in  order  to  develop  more  general  core  models.  In  particular,  a  core  model  which  predicts 
mass  transfer  rather  than  momentum  transfer  is  essential  for  the  characterization  of  diffusive  reacting  Hows 
such  as  chenhcal  lasers. 


C.  TWO-PHASE  MIXING  STUDY 

Powdered  fuels  hold  potential  advantages  for  advanced  propulsion  systems,  hut  their  successful 
utilization  requires  an  improved  understanding  of  the  complicated  injection,  mixing  and  combustion 
processes.  Some  of  the  fluid-dynamical  aspects  of  gas-particle  mixing  are  presented.  The  behavior  ot  a 
two-dimensional  jet  injected  into  .1  subsonic  cross  flow  was  studied  experimentally  tor  injection  ot  gases 
cf  different  densities. 

Solid  particles  may  require  a  relatively  long  residence  time  in  a  combustion  chamber  to  complete 
their  burning.  Injection  in  the  upstream  direction  appears  promising  because  of  increased  residence  time. 
For  a  constant  inlet  How  rate  application,  such  as  the  supersonic  ramjet  with  subsonic  cot  c  i.Uion.  the 
losses  corresponding  to  upstream  injection  are  exactly  compensated  by  reduced  shock  losses.  The  increased 
mixing  losses  caused  by  upstream  injection  were  also  investigated  for  flow  typical  of  scramjets  and 
stationary  power  plants. 

Particle  slip  relative  to  the  oxidizer  gases  was  demonstrated  to  increase  the  burning  rate.  The 
diffusion  limitation  which  often  controls  particle  combustion  may  be  reduced  if  sufficient  slip  can  be 
produced.  Combustion  instabilities  often  are  sustained  by  pressure  waves  in  the  system:  a  survey  of  the 
properties  of  various  waves  in  gas-particle  mixtures  was  prepared  for  publication. 

D.  MODELING  OF  MEAN  VALUE  DATA 

The  extension  of  universal  jet  width  growth  and  centerline  decay  laws  to  compressible  (variable 
density)  coaxial  jets  was  demonstrated  to  be  inadequate  for  characterization  of  this  class  of  How.  Data 
analyzed  showed  that  the  decay  exponent  of  the  centerline  mass  fraction  cannot  be  correlated  with 
either  the  ratio  of  the  initial  jet  velocity  to  that  of  the  free  stream,  or  the  ratio  of  the  initial  jet  mass  llux 
to  that  of  the  free  stream.  In  addition,  the  validity  of  the  Reichardt  inductive  theory  for  application  to 
compressible  flows  was  evaluated  and  found  to  predict  trer  Is  which  are  contrary  to  experimental  data 
when  a  Debt  gas  is  injected  into  a  slower  moving  heavy*  gas,  or  when  a  heavy  gas  is  injected  into  a  faster 
moving  light  gas. 
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II.  OBJECTIVE  AND  GOALS 


The  overall.  long-range  objective  of  the  Bel!  program  in  High  Speed  Mixing  and  Combustion  is 
to  provide  a  firm  technological  basis  for  the  design  of  high  energy  chemical  lasers,  and  advanced  propul¬ 
sion  systems,  using  both  gaseous  and  metal-loaded  fuels. 

The  goals  of  this  program  are: 

a.  To  develop  a  semiempirical  expression  for  the  Reynolds  momentum,  mass,  and  energy 
transport  applicable  to  each  of  the  flow  regions  of  reacting  and  nonreacting  turbulent  jets 
ti.e.,  from  the  injection  through  similarity  regions  I  and  for  ducted  flows.  This  modeling 
will  apply  in  the  core  region,  thereby  avoiding  *hc  usual  requirement  that  calculations  be 
started  downstream  of  the  core.  This  feature  of  the  modeling  is  important  since  often  rhe 
core  region  is  of  greatest  interest  in  hardware  design  consideration. 

b.  To  investigate  those  parameters  that  most  influence  the  dispersion  of  metal  particles  in  a 
non-tangential-injeclion  scheme.  An  optimum  particle  size  and  injection  mode  probably 
exists  for  which  the  time  from  injection  to  completion  of  combustion  is  a  minimum.  Slip 
between  particles  and  the  gas  tends  to  increase  the  burning  rate  by  elimination  of  the 
oxidizer  diffusion  limitation.  Methods  to  maintain  slip  in  a  combustor  by  suitable  modi fiea 
tion  of  the  flow  field  also  will  be  explored. 

e.  1  o  perform  an  experimental  and  theoretical  study  of  turbulent  mixing  downstream  of  a 
radial  injector  and  to  survey  experimental  and  analytical  mixing  literature  as  a  company- 
supported  effort. 
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III  PROGRAM  SUMMARY 


V  IURBULI  NT  MINIM'.  AND  RI  ACTING  I  LOW  (  II AKAC  1 1  RI/A I  ION 
I  INTROIIU TION 

TIk*  ability  to  nuke  meaningful  calculations  applicable  to  rocket  and  ramjet  combustors  and 
high  energy  lasers  has  been  a  long-range  goal  ot  designers  tor  a  number  of  years.  Unfortunately,  this  goal 
cannot  be  realized  until  prediction  ol  simultaneously  turbulent  mixing  and  reacting  flows  can  be  made 
with  confidence.  C  omputational  techniques  such  as  those  using  advanced  generalized  finite  element  met¬ 
hods.  e  g..  Rets.  I  and  2.  are  currently  available  tor  performing  such  calculations:  however,  the  lack  of 
adequate  understanding  ot  the  extremely  complex  physical  processes  and  their  interactions,  so  that  they 
can  be  modeled.  has  precluded  detailed  characterization  of  the  turbulent  flows  in  practical  systems.  That 
is.  before  n  eaningt ul  calculations  can  be  made,  the  turbulent  transfer  of  m  iss,  momentum,  and  energy, 
m  addition  to  chemical  reactions  kinetics,  and  their  complex  interactions  must  be  successfully  modeled. 
Flic  res*’  i:  ..1  task  reported  herein  has  the  ultimate  objective  of  improving  understanding  of  such  Hows  and 
develop  ine  appropriate  practical  models  and  computational  techniques. 


In  order  to  specify  turbulent  Hows.  Reynolds  transport  terms  which  characterize  the  turbu¬ 
lent  diffusion,  must  be  defined  for  each  of  the  transport  processes  in  order  to  solve  the  turbulent  species 
diffusion,  momentum,  and  energy  equa turns.  Fach  ot  these  processes  is  important  for  the  computation 
miicc  they  occur  simultaneously  and  influence  each  other:  however,  as  demonstrated  herein,  prediction 
ot  the  transfer  ot  mass  and  energy  are  more  critical  for  combustor  applications  than  the  prediction  of  the 
t runs t or  ot  momentum.  Apparently,  this  tact  frequently  has  been  overlooked  since  the  majority  of  tur¬ 
bulent  mixing  modeling  work  has  been  concerned  with  prediction  of  momentum  transport  alone,  e.g.. 

Ref.  .v  For  those  applications  in  which  the  transfer  of  mass  and  energy  are  considered,  the  assumption 
is  almost  always  made  that  turbulent  Schmidt,  Prandtl,  and  hence  Lewis  numbers  are  constant  ( frequently 
lunty ).  However,  a  number  ot  studies  has  shown  that  these  parameters  are  nut  constant  throughout  the 
mixing  region,  and  if  they  are  approximated  by  a  constant,  its  value  depends  at  least,  on  the  initial  con¬ 
ditions.  e.g  .  Scy  has  been  show  n  to  vary  several  fold.  Refs.  4  and  5. 


The  first  part  of  this  section  of  the  report  characterizes  the  relative  significance  of  the  various 
t>  pcs  <<  transport  in  turbulent  mixing  and  reacting  flows.  In  the  second  part,  turbulent  mass  l£)  and 
momentum  u  )  transfer  coefficients  are  presented  for  a  range  of  flow*  conditions,  these  coefficients  were 
obtained  by  direct  differentiation  of  experimental  data  using  a  procedure  designated  the  Inverse  Solution 
k.cciu;:v  I  he  third  portion  discus>es  the  extension  of  the  Inverse  Solution  Tv\  huiquc  to  chaij<.len/a'ion 
•>!  simnli  e  ussly  mixing  and  reacting  Hows,  and  the  List  poiuon  present'  an  appto.iji  lor  modeline 
practic.  i  mixing  and  reacting  Hows,  such  as  those  which  characterize  chemical  lasers. 


SIGNIFICANT  H  OF  Tl  RBI  Lb  NT  TRANSPORT  IN  MIXING  AND  Rh  AC  TIM;  I  LOWS 

Obviously,  m  a  combustor  or  a  cncmicai  laser  employing  diffusion  flames,  propellants  cannot 
react  before  they  are  mixed.  Ot  course,  their  residence  time  m  the  reaction  zone  is  an  important  factor 
in  the  determination  ot  the  completeness  of  t fie  reaction,  finally,  the  temperature  ot  the  reaction  zone 
also  is  important  since  chemical  reaction  rates  can  increase  exponentially  with  temperature  These  state¬ 
ments  appear  to  suggest  that  the  transfer  of  mass,  momentum,  and  energy  each  may  be  equally  important: 
however,  as  demonstrated  below,  such  a  coiidu\i*»n  definitely  is  not  warranted. 
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a.  C  omparison  of  Transport  Processes  for  C  oaxial  Following  Streams 

In  order  to  (|iijntitut'vcly  assess  the  relative  importance  ol  the  various  transport  pro¬ 
cesses.  a  series  of  computations  was  made  for  a  hydrogen  jet  mixing  and  reacting  with  air.  For  these 
computations,  an  explicit  finite-difference  numcncal  integration  technique  was  employed,  in  which  the 
governing  shear-layer  equations  were  solved  in  the  von  Miscs  coordinate  system.  The  finite-rate  kinetics 
were  included  in  the  calculations  using  a  last-running,  implicit,  quasilincari/cd  kinetic  subroutine.  A 
discussion  ot  the  computational  technique  is  presented  *n  Appendix  B.  For  these  calculations,  the  mithd 
velocity  and  hydrogen  concentration  profiles  were  assumed  to  be  identical  to  C  ase  I  A.  Chriss.  Ref.  6 
However,  in  order  to  obtain  reaction  rates  appropriate  for  combustors,  the  total  temperature  at  the  initial 
station  for  both  the  central  hydrogen  jet  and  the  air  stream  was  increased  from  C  hriss  *  reported  cold 
How  values  of  550  and  650°R.  respectively,  to  2025°R.  A  constant  eddy  viscosity,  e  -  0.02  Ibm/ft-scc 
was  used  in  the  calculation  as  a  representative  value,  together  with  turbulent  Schmidt  and  Prandtl  numbers 
of  0.0  (and  hence  a  Lewis  number  of  1 .0).  which  were  found  to  be  appropriate  in  Ref.  6.  Constant  values 
of  the  transport  coefficients  were  used  so  that  comparisons  of  the  various  predictions  would  be  independent 
of  the  shape  of  the  mean  profiles,  which  differ  greatly  for  the  various  cases.  If  an  eddy  viscosity  model 
that  was  dependent  on  the  local  flow  field  had  been  used,  values  of  e  would  have  differed  significantly 
in  each  of  the  four  cases.  By  using  constant  values  of  the  transport  coefficients,  the  relative  importance 
of  the  rate  ol  transport  of  mass,  momentum,  and  energy  were  dec  jpled  from  the  magnitude  and  shape 
ot  the  mean  profiles  making  direct  comparison  possible. 

Of  course,  using  this  procedure  would  not  be  expected  to  yield  results  consistent  with 
actual  experiments  because,  as  demonstrated  in  Appendix  A,  transport  coefficients  vary  throughout  the 
How  field  even  for  nonreactive  tlows.  Nevertheless,  quantitat;  e  determination  of  the  influence  on  the 
predicted  flow  field  of  a  change  in  one  transport  coefficient  could  bi  nost  easily  assessed  from  these 
calculations.  Results  are  presented  in  lugs.  I  to  3;  the  solid  curves  in  each  case  were  obtained  without  a 
reduction  in  any  ot  the  transport  coefficients.  As  indicated,  the  mass  (£j  and  momentum  le!  transport 
coefficients  were  each  decreased  in  turn  by  a  factor  of  10.  The  ratio  of  the  thermal  conductivity,  *,  to 
the  mean  specific  heat.  Cp,  also  was  varied  by  a  factor  of  10  since  this  ratio  is  dimensionally  consistent 
with  £  and  e.  i.e„  its  units  in  the  Hnglish  system,  are  Ibm/ft-sec.  Note,  thi.t  the  ratio  of  the  C’s  for 
K/Cp  -  0.0222  Ibm/ft-sec  and  Cp  for  k/c p  =  0.0022  Ibm/ft-sec  was  banded  between  0.77  and  L07;  the 
variation  in  cp  is  plotted  in  Fig.  *  for  the  cases  presented  in  Figs  1  to  3.  This  figure  shows  that  k  was 
^varied  by  at  most  10.7  and  at  least  by  7.7.  Therefore,  although  the  results  contained  some  influence  of 
c  variation:  this  influence  was  not  very  significant,  so  that  the  comparisons  in  Figs.  1  to  3  illustrate 
tne  significance  of  varying  both  k/c^  and  k. 

Calculations  presented  in  Figs.  I  to  3  were  made  starting  at  ?  =  0  with  initial  step 
(slug)  profiles  tor  V.  U.  and  Ty  using  the  procedure  described  in  Ref.  7.  Figure  l  shows  the  effect  ot 
these  variations  on  predicted  mass  fraction  water.  The  computed  water  profiles  at  various  axial  stations 
difter  gieatK  from  the  others  for  the  case  of  £  =  0  0022  Ibm  lt-sec.  hut  they  changed  very  little  forthe 
cm  responding  reductions  in  t  and  k  Cp.  Tliev  results  are  probably  the  most  important  of  all  the  com¬ 
parisons  because  they  clearly  demonstrate  that  the  extent  of  the  reaction  is  more  dependent  on  the  rate 
of  mass  transport  than  on  either  the  rate  of  momentum  or  energy  transport.  This  influence  is  not  con¬ 
fined  to  the  /  =  1  in.  station  but  persists  throughout  the  flow  field  up  to  /  =  7.3  in.,  the  final  axial  station 
considered. 


Axial  velocity,  U,  is  plotted  in  Fig.  2  for  the  same  set  of  input  parameters  used  in  Fie.  1 
This  figure  shows  that  both  e  and  £  influence  the  velocity  profiles  significantly  more  than  «/Cp.  at  least 
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Mgurc  3.  fcllect  of  Transport  Coefficients  on  the  Predicted  Static  Temperature  in  a  Reacting 
Hydrogen  Jet  at  i  =  1.0.  4. 1  5  and  7.3  in.;  (  ale illations  started  at  i  *  0.0 
U:  =  3300  It /see;  IL  =  525  fl/sec;  Tt  =  Tt  =  2025°  R 
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tor  /  <  4  in  In  Fig.  2.  the  e  =  0.002  Ihm  Tt-scc  case  exhibits  much  higher  centerline  velocities  because 
momentum  transport  is  inhibited  by  the  low  value  ol  e.  Interestingly,  the  £  0  0022  Ihm 'tt-scc  case 

caused  a  reduction  m  the  velocity  below  the  basis  case  because  the  density  gradients  are  much  steeper 
lor  this  case,  as  shown  m  Fig  5.  (Examination  of  the  momentum  equation  demonstrates  that  the  axial 
velocity  gradient  varies  inversely  with  the  local  density.) 

The  influence  of  the  various  transport  coefficients  on  static  temperature.  Ts.  is  presented 
in  Fig.  3.  Variations  in  £  and  */cp  are  shown  to  be  significant  in  the  determination  ol  the  Ts  profile 
Hardly  any  reaction  occurs  at  z  =  I .0  in.  for  the  £  =  0.0022  Ihm  tt-scc  case  since  without  mass  mixing, 
chemical  reaction  cannot  occur.  The  effect  oft  on  Ts  is  relatively  minor  as  was  the  case  for  the  water 
concentration  in  Fig.  I .  The  effect  of  decreasing  */cp  is  to  create  a  hot  spot  near  the  stoichiometric 
ignition  point. 


Figure  6  is  a  plot  summarizing  the  ignition  characteristics  for  each  of  th;  cases  presented 
in  Figs.  I  to  3.  The  ignition  characteristic  is  defined  as  the  maximum  value  of  total  temperature,  Tj  max 
minus  the  free-stream  total  temperature,  T-p*,.  attained  at  a  given  axial  station  in  the  flow.  Results  are 
very  similar  to  those  previously  presented.  A  10-fold  reduction  of  e  and  */cp  makes  relatively  little  dif¬ 
ference  in  Tj  max  except  to  shift  the  ignition  point  slightly.  However,  for  £  =  0.0022  Ibm/ft-sec,  Tj  niax 
is  reduced  by  almost  !000°K  (1800°R)  even  after  ignition  occurs,  and  the  ignition  delay  is  increased 
nearly  50'?.  Again  these  results  demonstrate  conclusively  the  fact  that  combustion  cannot  be  initiated 
until  significant  mass  mixing  has  occurred,  and  mass  mixing  depends  primarily  on  the  magnitude  ol  the 
mass  transfer  coefficient.  £ 

b.  Effect  of  C  oncentration,  Velocity,  and  Temperature  on  Chemical  Reaction 

A  second  series  of  calculations  was  made  to  demonstrate  the  importance  of  velocity, 
and  static  temperature  on  downstream  profiles  and  on  the  extent  of  the  H2/O2  chemical  reactions,  again 
using  the  computational  procedure  discussed  in  Appendix  B.  Comparisons  were  made  for  both  reacting 
and  nonreacting  cases.  The  limiting  cases  of  no  transport  of  either  mass,  momentum,  or  energy  were 
investigated. 


1 )  Mass  Transport 

No  additional  computation  was  required  to  demonstrate  that  when  no  mass  transport 
occurred,  no  chemical  reaction  could  occur.  Results  presented  in  Fig.  6  demonstrated  that  the  relative 
magnitude  of  £  compared  to  e  and  */cp  was  critical  from  the  standpoint  of  both  ignition  and  combustion. 
Obviously,  if  no  mixing  of  reactants  (mass  mixing)  occurs  at  all,  there  can  be  no  chemical  reaction  since 
constituents  must  be  intimately  mixed  at  the  molecular  level  before  chemical  reactions  may  occur  N<> 
mass  mixing  is  equivalent  to  each  stream  remaining  completely  unmixed,  and  therefore  unreacted.  Of 
course,  the  local  velocity,  temperature  and  initial  jet  and  free  stream  conditions  can  have  no  influence 
on  the  chemical  reaction  rate  until  mass  mixing  is  initiated.  The  case  of  no  mass  transport  at  all  is 
equivalent  to  £  =  0:  that  is.  initial  concentration  profiles  remain  unchanged  with  axial  position. 

2)  Momentum  Transport 

The  significance  of  momentum  transport  on  chemical  reaction  is  again  evaluated 
in  Fig.  7a  in  which  cases  with  and  without  momentum  transport  are  compared  In  one  case,  a  normal 
cold-flow  mixing  calculation  was  made  starting  with  experimental  hydrogen  concentration,  velocity. 
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Figure  6.  Fifed  of  Transport  Coefficients  on  Ignition  in  a  Heading  Hydrogen  Jet. 
Calculations  Started  at  i  =  0.0 
Uj  =  3300  ft /sec:  Ue  =  525  ft /sec;  Tj  =  Tjtf  =  2025°R 
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gun;  7a.  Importance  of  Momentum  Transfer  on  Predicted  Centerline  Hj  Mass  Fracti 
Ptedictions  Made  from  l  =  3.48  in.  using  e(r).  Le.-  =  1 .0  and  Scr  =  0.9.  I)  =  0.50  in. 


\ 


Bell  Aerospace  Company 


and  total  temperature  pioliles  from  z  =  3.48  m.  downstream  of  the  injection  station  for  C  ase  I  A.  of  <  Imss 
Kef.  (>.  Values  of  Scy  =  and  ley  =  1 .0  suggested  in  the  analysis  of  Kef  <>  were  used  in  the  calculations. 
The  transport  coefficients  used  in  the  computations  were  the  Composite  Trends  presented  in  the  next 
section  obtained  using  the  Inverse  Solution  Technique.  Experimental  centerline  data  points  are  plotted 
♦or  comparison.  In  the  second  ease,  the  velocity  was  not  permitted  to  vary  Irom  its  initial  valui,  i.c..  llu 
velocity  profile  was  "frozen"  at  its  value  at  3.48  in.;  therefore,  momentum  transport  was  completely 
precluded.  In  this  ease  (dashed  line  in  f  ig.  7a)  no  change  in  velocity  occurred  throughout  the  calculation; 
that  is.  the  velocity  throughout  the  flow  field  rein  lined  constant.  Results  show  that  very  little  difference 
occurred  in  predicted  values  ol  the  centerline  hydrogen  mass  fraction  for  these  two  cases,  both  of  which 
agreed  reasonably  well  with  the  experimental  data.  C  learly,  even  the  limiting  case  of  no  momentum 
transport  had  little  influence  on  the  predicted  concentration. 

Radial  experimental  profiles  for  this  ease  at  /. -  4.1$  and  7.30  in.  are  presented  in 
Figs.  7b  and  7c,  in  which  the  experimental  data  is  plotted  as  a  curve  and  the  computed  points  are  plotted 
as  symbols.  These  results  again  demonstrate  that  the  influence  of  momentum  transport  has  very  little 
influx  nee  on  predicted  concentration  profiles.  These  results  completely  confirm  those  of  the  last  section 
concerning  the  relative  importance  of  momentum  transfer. 

3)  Energy  Transport 

The  same  Chriss  I A  data  (Ref.  6)  used  in  the  last  Section  also  was  used  for  the 
evaluation  of  the  influence  of  initial  temperature  on  mixing  and  reacting  fWihr  Hows.  In  this  case,  step 
(or  slug;  velocity  and  concentration  profiles  were  used  to  initiate  the  calculations  from  the  injection 
station.  7=0.  so  that  the  ignii.  m  process  as  well  as  the  downstream  mixing  and  combustion  could  be 
modeled,  again  using  a  constant  value  of  e  =  0.02  Ibm/ft-sec  together  with  Sey  =  0.9  and  Ley  =  1.0. 
Chemical  reaction  rates  often  increase  exponentially  with  temperature;  therefore,  determination  of  the 
effect  of  a  1 80°  R  ( I00°K  >  variation  in  initial  temperature  from  2700  to  2880°  R  (1500  to  I600°K)  was 
of  interest.  Results  are  presented  in  Figs.  8  and  9.  Fig.  8  shows  that  the  minimum  difference  in  the  total 
temperature  rise  at  an  axial  station  0.2  in.  downstream  of  the  injector  was  only  about  90  R  (50  K).  Such 
a  relatively  small  difference  in  temperature  had  very  little  influence  on  the  mass  fraction  of  water  profiles 
shown  in  Fig.  0.  This  axial  station  was  selected  for  comparison  since  ignition  had  not  occurred  at  7=0. 1  in. 
for  both  cases,  and  the  Ty  and  water  concentration  differences  were  less  pronounced  at  z  >  0.3  in.  than 
at  7=0.2  in. 


Of  course,  further  reduction  in  the  initial  temperature  to  2025°  R  ( 1 125°  RE 
(Hi/Oi  stoichiometric  ignition  limit  is  830°K  at  1  atm),  did  cause  very  significant  effects  in  both  the 
total  temperature  rise  and  water  concentration.  In  this  case,  ignition  did  not  occur  until  about  7=0.7  in. 
These  results  are  presented  in  Figs  10  and  11.  Comparison  of  Figs.  8  to  1  1  suggest  that  at  least  for  the 
Hvair  diffusion  flame,  the  influence  of  temperature  on  ignition  delay  is  of  greater  significance  than  is  its 
influence  on  reaction  rates. 

3.  DIRECT  Dl .TERMINATION  OF  TURBULENT  TRANSPORT  COEFFICIENTS 

A  basic  understanding  of  turbulent  mixing  in  high  speed  flows  is  important  for  a  wide  range 
of  current  applications,  including  rocket  engines,  air-augmented  rockets,  supersonic  combustion  ramjets, 
and  chemical  lasers.  Prediction  of  jet  noise  also  requires  prediction  of  mean  and  turbulence  parameters 
in  coflowing  jets.  Unfortunately,  no  general  approach  exists  for  modeling  these  Hows  because  the  time- 
averaged  equations  of  change  cannot  be  solved  unless  experimental  information  is  available  concerning 
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Figure  7b.  Effect  of  Axial  Velocity  Variation  with  Distance  Downstream  on  Predicted  Hydrogen 
Mass  Fraction  Profile  at  z  =  4.1 5  in.  Calculations  Started  at  z  =  3.48  in.  using  Composite  Trend.  e(r). 

Case  1  A,  Chriss,  Ref.  6 
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Figure  7c.  F.ffect  of  Axial  Velocity  Variation  with  Distance  Downstream  on  Predicted  Hydrogen 
Mass  Fraction  Profile  at  z  =  7.30  in.  Calculations  Started  at  z  =  3.48  in.  using  Composite  Trend.  e(r) 

Case  1  A.  Chriss.  Ref.  6 
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Figure  1 0.  Effect  of  Initial  Temperature  on  Total  Temperature  Rise  for  a  Reactine  Hydrocen  Jet 
at  /  =  0.20  in:  Calculations  Started  at  /.  =  0.0  using  e  =  0  02  I  bm  ft -see.  LeT*  1 .0  and  Sct  =  O 
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,  !R  k  ,  !'rbu,c",ic-  a>  ,bm'  arc  a  wffidcnl  number  of  independent  equation*  to  uniquely 
‘V  Urh"  C.n'  ,ow'  rhuf  lv  experimental  data  required  lor  specification  of  the  Reynolds  trails-. 

'  r‘  ,i;r"'s-  J"d  ,or  1  "s  rcaso"  Prediction  ot  turbulent  flows,  by  whatever  approach  used,  must  be  con- 
mu  civil  scmicmpirK.il, 

.  ....  .  .  R1l‘1JSO",lblc  sl"-tcss  ,hl‘  Prediction  of  momentum  transfer  has  been  achieved  with  a  variety 

.  W\T  tt  d,y  V,S?M'y  J,'d  “,rbukncc  k,ntf,'c  energy  models,  e.g..  Kef  .  3  These  models  gwerallv 
sElnT".  ,RTdV  ‘‘T  ;,PP,K'a,'°"  tH  "",ll,sfK'cil*s  ‘’y  merely  assuming  that  the  turbulent  ' 

r  Jf!  T,bm  arc  CMns,ant-  S*'tb  a"  assumption  permits  computation  of  all  the 

Reynolds  transport  terms,  and  hence,  computation  of  the  turbulent  mass,  momentum,  and  energy  tran- 
sport  .  However,  as  discussed  earlier.  this  assumption  ot  constant  (and  generally  arbitrary!  ratios  between 

propa*a,,<’n  'naxs-  momentum. and  energy  (which  equals  the  ratios  of  the  appropriate  traa- 
sport  coefficients)  frequently  is  not  valid. 

...  t  n  !!  !*  un,lkcl>  *ha‘  a  “l,niversar  furbulence  models  applicable  to  all  types  of  practical  geometries 

for  wh  h"°  ,,0.nS  CVef  Ca"  **  dCVCk,ped  A,so’  limiiin*  thv  ran^‘  of  'l‘>w  conditions  aild^Lic  tries 

wh,^  modcl  ni,KUpply  |o  those  of  practical  interest  will  result  in  superior  predictions  within  this 

heTranlm.oZ' °  .  dCVC,°P  mixi"P  m,jJds  arplicable  to  specific  (lows  for  prediction  of  each  of 

the  transport  processes,  a  I*  chn.que  designated  the  Inverse  Solution  Technique  was  developed.  Ref  H  in 

an  I  Jh  oiTr?1  mean  concentration,  velocity,  and  temperature  profiles  are  differentiated  numerically. 

S  nC  tti,  annrn  .TVn  ,  eqUa,,°nS  °f chan*e  cvalua*td-  «cept  for  the  transport  coeffTcerrts. 

tnaLelv  niT  rnold?  frans,v,rt  terms  for  mass,  momentum,  and  energy  can  be  determined 

of  the  turb'd  n  cT' o  rd  ,°  C°"S,rU1Ct1  mixin*  modc,s  Wl,b<)ut  a,,y  assumption  concerning  the  nature 
ot  int  turbulent  Schmidt.  Lewis,  or  Prandti  numbers. 

minimi  h..  L^r  app,K'a'!T  of  ,hc  Solution  Technique  is  a  convenient  means  for  deter¬ 

mining  ihc  Reynolds  Transport  Terms.  Hot-wire  and  laser-Doppler  anemometry  techniquespem.it 

irect  measurement  ot  the  Reynolds  Transport  Terms  for  momentum  transfer  (the  Reynolds  Stresses  I  in 

2tim  rel  lurbu,encc  ,ntens,t,cs-  and  various  time  and  space  correlations,  and  laser-Roman  sattenng 

boi  I)1  rediou  s  a  nd't  im  '■*  corT  i'  *  d,f,CT,nina,i“  of  mass  transkr  efficients.  However,  such  measurements  are 
both  tedious  and  time-consuming  compared  to  pitot-static  probe,  gas  sampling,  and  stagnation  temperature 

Rev ToldTsT tS'  Wk'  h8rC  suffic*ent  tor  determining  each  of  the  transport  coefficients  (including  the 

w  s^tntll  mldsex  T 7 TV  ,  nVerse,So,"tiion  Technique.  A  detailed  discussion  of  this  technique 
„r ",  P  .d  R  K *  •,nd  9  Resu,,s  conclusively  demonstrated  that  meaningful  numerical  differentiation 
ol  data  can  he  accomphshed  if  the  data  are  properly  treated.  " 

Kef  10  for  ,^r  a"  aPPl,iCati0n  °,1  ‘he  '"VCrSC  So,M,ion  Technique,  mean-velocity  data  obtained  by  Moon, 
r  nH  ,h  1  aca,r  JCK  usm*  a  ho*  wire  anemometer,  were  used  to  determine  both  the  Reynolds  stress 
r.  and  the  eddy  viscosity  note,  that  r  =  e  iW  hr).  so  that  both  r  md  t  mas  he  considered  to  be 
momentum  trans.er  coefficients.  The  radial  r  profiles  are  plotted  he  I :  for  various  aSalsuLs 

J,Wf  W"  ’  tlK;  Kl'ynukls  s,n?'scs  determined  directly  from  the  x-w.re  measurements.  Agreement  is.  fuite 
».ood.  with  only  a  slight  displacement  of  the  maximum  values,  ds.nonstratme  that  valid  Experimental 
man,  velocity,  concentration,  and  static  temperature  profiles  may  be  used  in  a  similar  manner  for  direct 

d  wcZm  n."  f  t  f  T1<>,dS  "“"f  °ft  termS-  Such  a  determination  is  particularly  important  in  the 
development  of  models  lor  mass  and  energy  transport,  since  extensive  modeling  of  the  Reynolds  trans¬ 
port  terms  lor  these  processes  has  not  been  undertaken,  as  it  has  for  ease  for  momentum  tr^n 
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Dal;i  from  the  ten  other  cases  summarized  in  Table  I  were  analyzed  using  the  inverse  Solution 
luhnit|iic.  I  he  detailed  results  obtained  at  each  axial  station  lor  which  data  were  available*  together  with 
integral  injected  mass  and  momentum  balances,  are  presented  in  Appendix  A.  rigs.  A-l  to  A-IO.  (  areliil 
examination  of  these  results,  using  the  evaluation  procedure  developed  in  Kel  9,  indicated:  ! )  a  definite 
maximum  exists  tor  both  £  and  t  at  r  -'0,  2)  considerably  greater  variation  in  the  values  of  these 
cocMicionts  occurs  in  the  radial  (transverse)  direction  than  in  the  axial  (How)  direction.  3)  the  mass 
trail  si  ci  coefficients.  £.  (figures  designated  “a"),  exhibited  more  consistent  trends  than  do  the  momentum 
transfer  coclficients.  e.  (figures  designated  **b’);  therefore.  "Composite  Trends'’  were  obtained  from 
the  £  *  plots  rather  than  the  e-plots.  Of  course,  as  previously  demonstrated  £  results  also  are  of  greater 
interest  than  are  e  results.  These  simplified  C  omposite  Trends  arc  plotted  in  Fig.  13.  lor  the  data  sum¬ 
marized  in  Fables  I  and  2  in  which  the  transverse  or  radial  variation  of  the  mass  transport  coefficient. 

£,  was  retained,  while  using  an  average  axial  variation,  and  ignoring  those  regions  known  to  be  question¬ 
able  (see  Appendix  A).  Determination  of  trends  such  as  these  is  the  first  step  in  developing  mixing  models 
applicable  to  specific  (low  conditions  and  geometries.  This  technique  was  successfully  used  in  developing 
the  model  of  Ref.  4.  It  is  surprising  that  very  consistent  trends  were  obtained  for  such  a  wide  range  of 
experimental  data.  Fxcept  lor  C  urve  9,  Fig.  13  all  curves  have  very  nearly  the  same  shape  although  of 
course  their  amplitudes  depend  on  the  flow  conditions. 

Results  similar  to  those  obtained  for  the  mass  transfer  coefficients  summarized  in  Fig.  1 3 
were  obtained  for  the  eddy  viscosity,  e,  and  also  are  presented  in  detail  in  Appendix  A.  These  Composite 
e  *  plots  exhibited  the  same  trends  and  were  of  the  same  magnitude  as  were  the  £  -  plots.  Unfortunately, 
scatter  in  these  results,  caused  in  part  by  inconsistencies  of  up  to  207c  in  the  integral  balances  (see  Table  2). 
precluded  determination  of  valid  turbulent  Schmidt  numbers  (ScT  s  e/£).  Of  course,  when  valid  e  and  £ 
profiles  are  both  available,  their  ratio  at  a  point  is  equal  to  the  local  turbulent  Schmidt  number. 

Ihe  validity  ol  turbulent  transport  coclficients  can  best  be  assessed  by  using  them  to  predict 
the  llow  field  from  which  they  were  derived.  Ref.  8.  That  is.  the  coefficients  together  with  initial 
•clocitv.  concentration,  and  static-temperature  profiles  arc  used  in  a  numerical  integration  technique  to 
predict  downstream  profiles.  Agreement  of  these  predicted  profiles  with  the  original  experimental  data 
at  several  axial  stations  is  conclusive  evidence  that  the  transport  coefficients  are  consistent  with  these  data. 
Of  course,  such  agreement  does  not  assure  that  the  coefficients  are  “accurate."  They  can  never  be  more: 
accurate  than  the  original  experimental  data  from  which  the'  were  U.*rrvd;  however,  agreement  does 
guarantee  the  consistency  of  the  transport  coefficients  and  original  experimental  profiles. 

These  predictions  were  made  using  an  explicit  finite-difference  numerical  integration  pro¬ 
cedure  in  which  the  governing  shear  layer  equations  were  solved  in  the  von  Mises  coordinates,  see 
Appendix  B.  The  C  omposite  Trends  ol  Fig.  13  were  used  in  the  numerical  integration  in  each  case 
Since  momentum  transport  has  been  shown  to  be  less  significant  than  mass  transport,  and  in  some  cases 
appropriate  (  omposite  I  tends  for  e  could  not  he  obtained,  values  ol  e  equal  to  £  were  assumed  I  So r  -  |,i 
as  a  first  approximation  in  making  the  predictions  Ivpical  results  are  presented  in  Fig.  14  and  15  tor  the 
data  ol  (  ase  1 1  ,(  hnss.  Ret.  6,  and  C  ase  205KZawacki  and  Weinstein.  Ref.  I  I.  respectively.  Tabulated 
comparisons  for  these  cases  also  are  presented  in  Appendix  A.  Tables  A5  and  A7).  Agreement  between 
these  predictions  and  the  experimental  data  points  is  extremely  good  for  the  concentration  profiles  and 
reasonable  lor  the  velocity  profiles  as  well.  These  results  demonstrate  conclusively  that  the  Composite 
Trends  adequately  represented  both  £  and  e  for  these  cases.  Although  the  agreement  was  not  as  good 
in  every  case  presented  in  Appendix  A.  reasonable  results  were  attained  in  all  cases,  which  tend  to 
validate  the  Composite  Trends  of  Fig.  1 3.  Note,  starting  with  a  composite  trend  or  a  model  for  £  and 
assuming  that  e  *  £  is  not  the  same  as  predicting  e  using  the  turbulence  kinetic  energy  approach,  or  art  eddy 
viscosity  model,  and  assuming  £  =  e  as  is  frequently  done.  e  g..  Ref.  3.  although  Scj  is  unity  in  both  cases. 
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List  of  Symbols 


D;  Injected  Jet  Diameter 


M  -  Mach  Number 
0  -  Density 


U  •  Axial  Velocity 

(  >e  -  External 

(  )-  .  Injected 


Momentum  integral  balances  were  within  4%  of  their  average  value  at  any  axial  station, 

*  Injected  mass  Integra'  balances  v*re  within  20%  o#  the»r  average  value  at  any  axial  station 
and  v !*re  within  VJ%  of  their  average  value  at  80%  of  the  axial  stations. 
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COAXIAI  Jl  T  MIXINi:  DATA  ANALYSIS 
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Ligure  13.  Composite  Turbulent  Mass  Transfer  Coefficient 
Profiles  for  Data  of  Tables  1  and  2 
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Predictions  Made  Using  Composite  Trend.  £  (r).  and  Scy  =  Ley  =  1.0 
Data  of  Chriss.  Ref.  6. Case  IE.  Calculations  started  at  z  =  2.74  in. 
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Figure  15.  Predicted  and  Experimental  Velocity  and  Freon  Mass  Fraction  Profiles. 
Predictions  Made  Using  C  omposite  Trend.  £  (r),  Scy  =  Ley  =  1.0 
Data  of  Zawacki  and  Weinstein.  Ref  I  I ,  Case  205 F 
Calculations  started  at  z  =  2.0  in. 


Bell  Aerospace  Company 


4  TK  HNKU1  FOR  I  VALUATION  OF  MI\IN<;  MODI  I  S 

Previous  work  on  this  contract  resulted  in  the  development  of  scmtcmpirical  eddy  viscomix - 
t\  pe  mixing  models,  applicable  in  the  transition  region  of  circular  free  and  coaxial  jets  for  a  reasonably  wtue 
range  ol  flow  conditions,  Reis.  7,  13.  Ot  course,  neither  these  models,  nor  most  likely  any  other  semiempu- 
ical  mixing  model,  ever  will  be  applicable  to  ail  turbulent  jet  flows.  Therefore,  determination  of  the  range  ol 
conditions  for  which  a  particular  model  is  applicable  is  extremely  important.  The  Inverse  Solution  Techni¬ 
que  w  as  used  to  evaluate  the  applicability  of  the  models  of  Refs.  7  and  I  3  for  several  specific  Hows.  Clear  l\ . 
modification  ot  any  mode!  to  make  it  applicable  to  a  particular  type  of  flow  could  be  accomplished  using  the 
same  approach. 

For  example,  the  maximum  and  minimum  values  of  £  (obtained  at  the  initial  and  final  stations) 
for  Case  IF.Chriss  Ref.  6.  Case  205F,  Zawacki  and  Weinstein,  Ref.  1 1 ,  and  Case  I  A,  Chriss,  Ref.  6.  computed 
with  the  eddy  viscosity  model  ot  Ref.  7.  using  a  value  of  Ley  =  1 .0  and  the  values  of  Scy  indicated  in  each 
figure,  are  plotted  in  Figs.  16  to  18.  respectively.  The  axial  stations  at  which  the  predictions  were  made  are 
designated  on  the  appropriate  curves.  The  appropriate  Composite  Trend  of  Fig.  1 3  is  replotted  in  each  figure 
for  comparison. 

Agreement  between  the  values  of  £  obtained  using  the  model  of  Ref.  7  and  the  Composite 
Trend  is  not  very  satisfactory  in  Fig.  16.  Since  the  predictions  in  Fig.  14  using  this  Composite  Trend  closely 
predicted  the  experimental  data,  the  mode!  is  demonstrated  to  be  somewhat  inadequate  for  the  Chriss  11: 
data.  That  is.  the  magnitude  of  £  predicted  by  the  model  is  too  low.  which  will  lead  to  underestimation  of 
rhe  mixing.  Fairly  close  agreement  between  the  Composite  Trend  and  the  model  of  Ref.  7  was  achieved  with 
the  Zawacki  and  Weinstein  data  in  Fig.  17.  Since  good  agreement  was  obtained  with  the  experimental  data  in 
Fig.  I  5.  using  the  Composite  Trend,  the  model  should  be  adequate  for  this  case  and  should  result  in  better 
predictions  than  for  the  previous  case. 

Comparisons  in  Fig.  18  illustrate  t he  use  ol  the  Composite  Trends  to  determine  which  of  two 
models  is  most  appropriate  tor  Case  I  A.  Chriss.  Ref.  6.  The  model  used  to  generate  curve  (a)  of  Fig.  18  was 
identical  to  that  used  to  generate  curves  (b).  except  for  a  factor  D  s  Ipu/pinjnl  ^  as  discussed  in  Refs.  7 

and  1 3.  Clearly,  improved  predictions  are  anticipated  when  including  the  factor  D  This  result  is  confirmed 

in  Fig.  19  in  which  Y  and  l1  profiles  at  z  =  7.30  in.  are  compared  with  and  without  the  D  The  agreement 
between  predictions  made  using  either  model  and  the  experimental  data  is  not  as  good  as  that  obtained  with 
the  C  omposite  Trend,  the  dashed  curve  in  Fig.  19.  This  result  w  as  expected,  since  the  Composite  Trend  was 
obtained  using  the  Inverse  Solution  Technique;  it  is  specific  to  a  particular  geometry  and  set  of  flow  conditions. 
A  mixing  model  attempts  to  correlate  a  fairly  wide  range  of  How  conditions,  and  therefore,  generally  will  not 
be  precise  for  a  particular  flow.  Clearly,  when  used  as  illustrated,  the  Composite  Trends  obtained  from  the 
Inverse  Solution  Technique  are  very  useful  in  the  synthesis  of  turbulent  mixing  models. 
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parison  of  Composite  Trend,  £  (r),  and  Predicted  $  Using  Model  of  Ref.  7  and  S 
LeT  =  1.0,  Case  205  Fr  Zawacki  and  Weinstein,  Ref.  1 1 


Predicted  (a)  Without  0_(o)  With  0 


l  ase  I A  (  hriss,  Kcl  . 


U/U-.  U.  =  525  ft/sec 
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u 


*  Predicted  with  Composite  Trend 
EupChnss  1A 

-  Predicted  Without 

-  Predicted  With 


r/a.  a  *  0.25  in. 


0  0.5  1.0  1.5  2.0  2.5  3.0 

r/a 

(b)  Mass  Fraction 


Pigurj  PC  id  Op  on  Predicted  Mean  Velocity  and  Ihdrneen  \Ciss  I ' rac t u > n  Profile 

at  i-  7.30  Inches.  Scj  =  0.9,  Lej  =  1.0.  Data  ol  Case  lA,C'hriss.  Ret.  6 
Calculations  started  at  z  =  3.48  in. 
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5  ANAI  YSIS  OF  MIXINt;  AND  REACTING  FLOWS 

woul  ii,.  ,  A,'rlK'"'0?  l,U>  '"T*  S<’,U"“n  Tcchm<|U0  “»  «mullanwusly  mixing  and  readme  Hows 
P1'  K  Realise  «»(  the  extreme  difficulty  ...  obtaining  duvet  measurement  uf.heKeynohK 

I  ransport  Terms  lor  mass,  mumcnfim.  and  energy  in  such  Hows  Ol  course,  ultimately  laser*, ppter  amt 
laser  Raman  scattering  teelim.|ues  may  make  such  measurements  feasible .  however,  gas  samntoig.  nitot- 
probe.  and  stagnation  temperature  measurements  are  tar  simpler  to  obtain.  Therefore,  for  the  foreseeable 
future  at  least.  .1  would  be  desirable  to  determine  Reynolds  Transport  Terms  directly  from  such  probe 
data  in  a  manner  similar  to  that  discussed  in  the  last  section  for  nonreactive  flows.  As  demonstrated  below 

''*rt  'I*  T  V"1’  '"TT  S°!U"0n  Tcc  hnk,UC-  at  leas‘  for  suth  as  «he  core  and  tram,,, on 

f  °*  '•1,  a",u,1larK,v  wh,ch  onl>  nu\wg.  and  hence  chemical  reaction  hasoccunwL.  AppTi- 

eation  of  the  technique  to  flows  which  have  reacted  significantly  has  not  yet  been  attempted. 

The  hydrogen-air  combustion  data  of  Cohen  and  Guile.  Ref.  14.  was  used  to  demonstrate  this 
procedure.  Table  3  is  a  summary  o»  the  integral  hydrogen  mass  and  momentum  balances  at  each  axial 

v  .rv'm ‘°r.W  n  ?PCrimCn,i"  da,a  WCfe  oh,ained-  B<,,h  ,he  hydrogen  flow  rate  and  momentum 

arv  considerably  Irom  station  to  station:  however,  the  variation  of  the  integral  balances  at  the  7  10 

and  14-mch  stations  is  minimal,  (although  the  injected  mass  balance  did  vary  from  +1 7  3  to  -7 Of  there- 
lore.  these  stations  were  used  exclusively  in  the  analysis.  In  these  experiments  in  which  only  limited 
quantities  ot  water  were '  formed .  detailed  concentration  velocity,  and  total  temperature  data. (obtained 
with  conventional  sampling  probes!  were  available.  The  data  for  the  four  axial  stations  considered  for  the 
ugh-temperature  vitiated-air  case  are  presented  in  Figs.  20a  to  20d.  respectively.  The  dashed  line  in 
Fig.  .0a  is  a  smoothed  curve  drawn  through  the  “actual"  experimental  data  at  z  »  14  in.  The«#A  curves 
were  obtained  by  adding  the  injected  hydrogen  consumed  during  combustion  to  the  measured  (actual! 
hydrogen;  the  technique  tor  computing  this  “total"  hydrogen  is  discussed  below.  Even  at  the  furthest 
station  downstream  </  =  14  in. )  at  which  the  maximum  combustion  occurred,  very  little  difference 
exists  between  the  measured  and  "total"  hydrogen  concentration  profiles:  therefore,  measured  hydrogen 
as  not  plotted  at  the  earlier  stations.  Note  however,  that  a  I'?  variation  in  hvdrogen  mass  fraction  due 
to  reaction  corresponds  to  a  variation  in  water  mass  fraction. 

In  the  analysis,  the  assumption  was  made  that  hydrogen  and  water  diffuse  at  the  same  rate 

!£;  ,a  r  ^e  sa,nc  Va,w  °*  l  "  appl,cable  ho,h  spe“ex  they  mix  with  the  remaining  constituents,  so 
hur  mass  mixing  can  be  characterized  Irom  the  total  hydrogen  profiles  alone.  However  the  chance 
m  velocity. 'temperature,  and  density  caused  by  combustion  due  to  the  heat  release  is  included  by  using 
I  tese  experimentally  determined  profiles.  Thai  is  the  assumption  is  made  that  the  water  formed  remains 
idiltusesi  with  the  hydrogen  as  though  it  had  not  reacted,  i.e  .  the  mass  mixing  is  characterized bv  “total" 
hydrogen.  Clearly,  this  assumption  is  valid  when  only  limited  reaction  occurs  since  in  this  case,  hvdroeen 
depletion  and  the  water  lonnahon  is  negligible:  its  jpP|,eah,li.x  ..occ  Ul!h  ,xu.„slve  reac,jon  must 

demonstrated  Note,  however.  that  H  i  i-  even  approximate!.,  the  -ame  tor  hvdrogen  and  water. the 
exhtmiue  could  be  applied  to  highly  reactive  flows.  Currently .  the  assumption  that  {  is  identical  for  all 
species  in  multispecies  turbulent  flows  is  almost  always  made.  e  g.  Ref.  15;  however,  its  validity  must  be 
determined  by  careful  analysis  of  detailed  data. 

••  -onvnrr .  f  !n  !he,prcscn' v  al- tsis-  ,h°  Watcr  ,hat  ,brmed  duc  to  combustion  of  the  injected  hydrogen  was 
converted  to  hydrogen,  h  u,  lr  .  using  the  relation 
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TABU  3 

CONSISTENCY  Of  REACT!  NO  COAXIAL 
HYDRCXiEN  JI  r  DATA* 


A*ial  Distance 
Downstream 
Inches 

Hydrogen  Flow 
Rate 

Ibm/sec 

%  Variation 
From  Mean 

Momentum 

Ibf 

%  Variation 
From  Mean 

0  085 

0.0310 

-54.C 

103 

21  6 

2.0 

0.0788 

+  17.3 

143 

+  7.5 

4  0 

0.0922 

+372 

134 

+  0.8 

7.0 

0.0711 

+  5.8 

136 

*  1-5 

100 

00675 

+  0.5 

146 

4  9.8 

14.0 

0.0625 

-7.0 

137 

+  3.0 

*Data  of  Cohen  and  Guile,  High  Air  Case,  Ref.  14 


Total  Hydrogen  Flow  Rate  *  fpUYyAI- 
Momentum  *  £  (pUUA). 

where  Yt  =  |YH  |  of  Eq.  (41 
^  Total 

and  t  designates  each  sampling  point 
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ul  Velocity 


High  Temperature  Vitiated  Air  Case 


f  igure  20d.  Experimental  Total  Temperature  Pro! lies  at  /  -  2.0,  7.0,  10.0  and  14.0  Incites 
Solid  C  urves  represent  Smoothed  Dated  of  Cohen  and  Guile,  Ref.  14. 

High  Temperature  Vitiated  Air  Case 
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Of  course,  the  validity  of  such  derived  transport  coefficients  can  only  he  demonstrated  hy 
using  (hem  in  the  numerical  integration  procedure  discussed  in  the  last  section,  in  which  experimental 
profiles  are  used  as  the  initial  or  starting  conditions  in  the  numerical  integration  of  the  appropriate 
governing  equations  (See  Appendix  B).  Predictions  made  using  this  technique  with  the  Composite  f  of 
Fig  21  are  compared  with  the  experimental  hydrogen  concentration  m  Pig.  22a,  the  axial  velocity  in 
Fig.  22h,  the  total  temperature  in  Pig.  22c  and  the  water  concentration  in  Fig.  22d.  Predictions  were  made 
with  the  assumption  of  both  frozen  flow  (no  reaction )  and  finite  rate  \\2  /()2  chemical  kinetics.  The  finite 
rate  kinetic  computational  technique  also  is  discussed  in  Appendix  B.  Predictions  made  using  both  the 
finite  rate  kinetic's  and  the  frozen  flow  assumptions  agreed  to  within  ♦  10%  of  the  experimental  hydrogen 
concentrations  (Fig.  22a)  except  at  the  z  ~  4  in.  station  for  r  <  0.4  in.,  in  which  the  experimental  hydro¬ 
gen  concentration  exhibited  what  appears  to  be  an  anomalously  steep  gradient.  The  agreement  of  the  two 
sets  of  predictions  is  not  surprising  since  relatively  little  hydrogen  was  reacted,  (See  Fig.  20a).  In  Fig.  22b, 
predictions  of  the  velocity  are  within  ♦  20%  data.  The  fact  that  the  predicted  velocity  is  in  poorer  agree¬ 
ment  than  the  concentration  is  not  surprising,  since  e  and  «/cp  were  assumed  equal  to  {.  Nevertheless,  the 
agreement  achieved  appears  very  adequate  for  design  and  evaluation  of  most  hardware.  Figure  22c  shows 
that  the  finite  rate  kinetic  calculations  predict  the  total  temperature  to  within  about  ±10%  of  the  experi¬ 
mental  value.  Obviously,  the  total  temperature  predictions  will  be  higher  for  the  reactive  case  (finite  rate 
kinetics)  than  for  the  nonreactive  case  (frozen  flow).  Nevertheless,  a  temperature  gradient  is  computed 
for  the  case  of  frozen  /low  because  the  free  stream  is  considerably  hotter  than  the  hydrogen  jet. 

The  predicted  water  profiles  in  Fig.  22d  do  not  agree  as  closely  with  the  experimental  data 
as  did  the  other  parameters  of  Fig.  22.  In  part,  the  discrepancy  may  stem  from  the  experimental  diffi¬ 
culties  of  quantitative  water  determinations.  However,  the  significant  and  consistently  higher  prediction 
of  water  obtained  using  the  finite  rate  kinetics  suggests  that  other  effects  also  may  be  partially  responsible 
for  these  discrepancies.  The  assumption  that  nitrogen  and  water  both  diffuse  at  the  same  rate  as  hydrogen 
in  Eqs.  (2)  and  (4),  so  that  their  mass  transfer  coefficients.  |j,  are  ail  identical,  must  be  questioned.  These 
are  the  most  likely  explanations  for  the  fact  that  the  Composite  $  =  e  =  k  /cp  used  in  the  computations  re¬ 
sulted  in  rather  good  predictions  of  each  of  the  transport  processes,  as  demonstrated  in  Figs.  22a.  h.  and  c. 
and  of  the  kinetics  as  well,  at  least  as  characterized  by  the  heat  release  in  Fig.  22c,  but  did  not  result  in 
as  good  predictions  of  the  water  profile  in  Fig.  22d.  The  possibility  is  discussed  in  the  next  section  that 
different  mass  transfer  coefficients  must  be  defined  for  at  least  some  of  the  important  reacting  species, 
in  order  to  adequately  characterize  turbulent  reacting  flows.  The  results  of  Fig.  22  tend  to  substantiate 
this  suggestion.  For  example,  if  O  is  assumccl  to  be  less  than  ,  an  overall  reduction  in  the  predicted 
water  level  would  result:  additional  data  arc  required  before  definite  conclusions  can  be  drawn. 
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Figure  21.  Experimental  am!  Smoothed  £  (r).  Mass  Transfer  Coefficient.  In  The  Transition 
Region  of  a  Reacting  Hydrogen/Air  Jet  (Data  of  Cohen  and  Guile.  Ref.  14 
High  Temperature  Vitiated  Air  Case). 
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Predicted  (Frozen  Flow)  Using 
Composite  £  (ri 


Predicted  (Finite  Rate  Kinetics) 
Using  Composite  £  (r). 

O  Exp. 


4.0  5.0  6.0  7.0  8.0 


U  •  Axial  Velocity  (Thousands  of  feet) 


Figure  22b.  Predicted  and  Experimental  Velocity  Profiles  Predictions  Made  from  /  =  2.0  In.  using 
Composite  £  (r),  Ley  =  1.0  and  Scy  =  1.0  (Data  of  Cohen  and  Guile,  Reference,  High 
Temperature  Vitiated  Air  Case). 
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O  Exp.,  Cohen  and  Guile 
*“*•**■•  Predicted  (Frozen  Plow! 

Using  Composite  £  (r) 

Predicted  (Finite  Rate  Kinetics) 
Us.ng  Composite  £  (r) 


Tj  Total  Temperature 


re  22c.  Predicted  and  Experimental  Total  Temperature  Profiles.  Predictions  Made  From 
Composite  £  (r).  Ley  =  1 .0  and  Scy  =  I  0  (Data  of  Cohen  and  Guile.  Ref.  14. 

Hiph  Temperature  Vitiated  Air  Case!. 
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H20  Mass  Fraction 
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ST  '  V,|;  is  "k>  water  that  was  formed  when  the  injected  hydro** t  ‘•burned  “ 

— -i  vh,  itrsiT*  . . -  -  - <•-—* "» 


iyh,oIv  *  ivHl0r 


rN, 

IVN,I. 


Cl 


|V"-0,„  •  IV».OlTola,  -  IVHl0l 


f3) 


SS^ftgSuC  <VH. » PSKIhe  too- 


IV. 


Hl  To'al  =  ”j  +  ,YH,Ic 
Tl«e  solid  curves  in  Fig.  20a  were  obtained  from  F„.  ,4 Density  profiles 

wore  generated  using  the  .elation 


f4) 


requited  for  the  calculations. 


PMW 

RTt. 


(5) 


«— 

second -order,  running  smoothiim  techniiiuc*  R'fs  R  ™  \  <>  oh  CL  mplished  u*inga  tive-ptunt. 

such  as  central  differences  r£S'  a^d  famH  ,  i  nTencal  J'^re„«ation  procedures 

hue  drawn  m  the  figure  represents  the  average  radial  vanafion  or  (  on  P,,  ,e  Trend  Sr t he^rV as  , 
stations  indicated.  Results  obtained  at  the  initial  station  /  =  Y  „JL“^  ,  a  u  h  ‘ 

a  large  number  of  neaative  values  oft  This  result  is  n  n  "■  1  1  1  Plotted  because  rhey  included 

obtained  at  intermediate  stations  are  valid  'he  ^  P°,ntS‘  S°  thaf 
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(>.  MODI  !  INC.  OK  MIXING  AM)  REACTING  FLOWS 

Analysis  of*  propulsion  and  chemical  laser  systems  require  the  ability  to  model  s  nmltaneousK 
mixing  and  reacting  turbulent  tlows.  Snell  analyses  have  been  ot  interest  tor  many  years  in  the  chemical 
process  industry  as  well,  and  various  approaches  have  been  attempted  A  recent  discussion  ot  some  of 
the  previous  wotk.  jIoih*  with  a  suggested  new  approach*  was  presented  In  Spalding  in  Ref.  lb  In  this 
new  approach,  a  single  simple  global  local  volumetric  reaction  rate  was  vlefined.  which  was  assumed  to  he 
limited  by  I  ho  rate  ot  break-up  of  lumps  (or  eddies!  of  unbumed  gases.  Spalding  stated  that  the  local 
turbulence  level  has  a  strong  influence  on  the  local  reaction  rates  and  that  in  most  eases  it  is  a  dominant 
influence.  This  suggestion  is  in  agreement  with  calculation  by  Donaldson.  Kef.  I  7,  who  demonstrated 
that  the  turbulence  level  can  exert  a  very  strong  influence  on  the  rate  of  bimolecular  chemical  reactions 
in  a  simple  idealized  system. 

This  influence  is  easy  to  explain  physically  by  noting  that  chemical  kinetic  reaction  rate  con¬ 
stants  are  applicable  to  instantaneous  rather  than  average  or  mean  concentrations  of  the  reacting  species, 
since  they  normally  arc  obtained  in  reacting  systems  at  rest:  therefore,  the  level  of  turbulence  which  in¬ 
fluences  instantaneous  concentrations,  must  influence  the  rates  of  reaction  as  well.  In  addition,  the  com¬ 
bustion  process  must  also  exert  an  influence  on  the  turbulent  mixing  because  it  creates  large  local  con¬ 
centration.  velocity,  temperature,  and  density  gradients,  which  directly  modify  the  local  turbulence  level. 

One  obvious  limitation  of  the  analysis  presented  by  Spalding  in  Ref.  16  is  the  fact  that  a 
single,  global  reaction,  characterized  by  a  “maximum"  value  of  the  volumetric  burning  rate  was  employed. 
This  value  depends  on  the  local  turbulence  level  (i.e..  How  field)  as  well  as  on  the  particular  set  of  re¬ 
actant's  employed.  The  ability  of  such  a  formulation  to  yield  quantitative  predictions  of  a  complex  sys¬ 
tem  consisting  of  dozens  of  reactions,  many  involving  excited  species,  such  as  are  required  in  chemical 
laser  analysis,  e.g..  Ref.  18.  must  be  questioned.  For  example,  it  is  very  unlikely  that  a  different  “maxi¬ 
mum"  volumetric  burning  rate  (which  depends  on  the  local  turbulence)  could  be  defined  for  each  of  the 
many  of  kinetic  reaction  steps  which  characterize  a  chemical  laser. 

Since  the  interaction  of  turbulent  mixingon  the  rate  of  chemical  reactions  and  the  influence 
of  the  reactions  on  the  mixing  is  extremely  complex,  it  is  worthw  hile  to  discuss  a  direct  approach,  which 
demonstrates  the  various  interactions.  For  the  bimolecular  reaction 


a  +  b  =  e  +  d 


the  forward  reaction  rate  may  he  written  as 


Forward  Reaction  Rate 


ui 


=  ^  Ya  Y„ 


where  Ku  is  the  usual  reaction  rate  constant,  and  Y  is  the  instantaneous  local  mass  fraction.  Since  these 
mass  fractions  are  fluctuating  quantities,  this  reaction  rate  is  an  instantaneous,  fluctuating  quantity  as 
well. 


Conceptually,  these  instantaneous  concentrations  could  be  measured  at  a  point  in  the  flow 
as  functions  of  time,  perhaps  using  lascr-Raman  scattering  techniques,  and  the  corresponding  reaction 
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mIcn  tide rm mod  directly  trom  Iq.  (61.  rime-averaging  over  an  appropriate  interval,  r,  would  vield  a 
moan  rate  ol  the  forward  reaction  applicable  to  that  point  in  the  How  (  Icarly.  this  mean  reaction  rate 
would  depend  on  the  turbulence  level,  since  Ihe  instantaneous  values  of  the  fluctuating  m**  tract  ion 
ol  the  reactants  are  directly  influenced  by  the  turbulence  as  well  as  the  chemistry'.  If  the  time-averaged 
concentrations  are  determined  in  the  conventional  manner,  an  average  forward  reaction  rale  “const airt  “ 
Ka  may  be  defined  by  Ihe  relation. 


Mean  forward  Reaction  Rate  ^ 


j 


(Forward  Reaction  Raleldt  ~  Ka  Ya 


(n 


Ot  course  the  term  Ka  will  no  longer  be  a  “constant**  (primarily  a  function  of  temperature);  as  is  fC3„  be¬ 
cause  it  is  dependent  on  the  turbulence  level  of  the  flow  field*  as  well  as  the  chemical  kinetics:  newt  he- 
less,  this  value  when  used  with  time-averaged  concentrations,  would  yield  the  precise  rate  of  reaction  ap¬ 
propriate  to  the  sampling  point  and  time  interval,  rt  selected.  This  discussion  demonstrates  the  com¬ 
plexity  of  the  problems  that  must  be  solved  when  considering  the  interrelation  of  turbulence  and  chemical 
reaction. 


Reacting  flows  such  as  chemical  lasers,  that  require  a  large  number  of  chemical  kinetic  re¬ 
actions  to  adequately  characterize  the  chemistry  must  utilize  existing  chemical  kinetics,  at  least  for  the. 
foreseeable  future,  i.e.  Ka  =  Ka,  which  simply  means  that  accommodation  of  the  interaction  of  turbulence 
and  kinetics  must  be  empirically  included  in  the  computation  of  the  flow  fieid,  rather  than  bv  reformu¬ 
lation  of  the  kinetics.  Therefore,  the  mass  mixing  coefficients  rather  than  the  Kasmu»t  be  modi¬ 
fied  so  that  the  time-averaged  concentrations  ol  reactants  computed  from  solution  of  the  time-averaeed 
form  of  the  species  diffusion  equations  will  yield  realistic  reaction  rates  when  used  with  the  conventional 
kinetic  rate  constants,  at  least  tor  the  critical  constituents  such  as  the  activated  species  in  a  chemical 
laser  cavity.  That  is.  the  Ya  s  and  s  computed  in  the  mixing  analysis  together  with  the  conventional 
Ka's  will  yield  the  appropriate  mean  forward  reaction  rate.  Of  course,  an  identical  analysis  applies  to  the 
reverse  reactions  as  well. 

The  feasibility  of  applying  this  approach  was  evaluated  using  the  Cohen  data  presented  in 
the  last  s  ection.  An  effort  was  made  to  modify  the  eddy  viscosity  model  of  Ref.  7  so  that  the  standard. 

H2 102  kinetics,  discussed  in  Appendix  B.  could  be  applied.  The  simplest  assumptions  were  made  for 
this  modeling:  ( I )  A  single  value  of  £  applicable  to  all  species  may  be  employed.  (2)  The  turbulent 
Schmidt  and  Lewis  (and  hence  Prandtl)  numbers  were  assumed  to  be  unity  so  that  =  (*!cplm  -^rn 
could  be  predicted  directly  from  the  model:  the  subscript  in  denotes  mixing  alone,  i.e.,  no  combust imt. 

An  empirical  factor.  C,  was  arbitrarily  defined  to  account  for  the  interaction  of  mixing  and  combus¬ 
tion  so  that 


*m+c  ~  rCp)m+c  =  tm+c  ^lm  <*► 

where  m+c  denotes  simultaneous  mixing  and  combustion  This  particular  form  of  Eq.  (81  was  selected 
because  of  Ihe  considerable  success  obtained  using  a  similar  form  for  the  correlation  of  practical  combus¬ 
tor  performance.  Ref.  19.  Of  course,  such  a  simple  approach  will  only  be  useful  if  it  can  be  related  to  - 
reacting,  turbulent  flows  in  a  relatively  simple  and  consistent  manner.  It  must  predict  the  charac¬ 
teristics  that  define  the  flow  such  as  Y,  U,  p.  Ty.  which  means  the  chemistry  and  heat  release  also  must 
be  adequate  described.  That  i>  resulting  values  of  |m+(.  must  vield  realistic  computation  oftocal  reaUant 
concentrations. 
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Computations  made  assuming  C  to  be  a  constant  equal  to  unity  resulted  in  prediction  ot 
significant  overmixing  of  both  hvdropen  concentration  and  velocity.  Comparison  of  values  of  $  obtained 
with  the  model  diowed  it  to  be  about  four  times  the  Composite  Trend  of  Fig.  -  I  Therefore,  a  value  of 
C  =  0.25  was  tried.  Results  of  the  computation  are  presented  in  Ftps.  23  and  24.  Comparison  of  the  re¬ 
sults  of  Fiji.  21 .  obtained  using  the  Inverse  Solution  Technique,  which  arc  rcplottcd  in  Fig.  23,  with  values 
°f  *m+c  r|L'J>c,ed  using  I  q  iM  and  (  *  0  25.  show  that  moral*  rather  good  agreement  was  achieved.  The 
predictions  for  the  axial  stations  /  =  7  and  10  in.  completely  bracket  the  Inverse  Solution  Technique  re¬ 
sults.  The  fact  that  the  prediction  ot  f  at  /  -  14  in.  is  high,  suggests  that  the  e-model  of  Ref.  7  contains 
more  axial  variation  than  needed.  e*en  when  this  variation  is  reduced  four  fold  by  assuming  C  =  0.25. 

The  comparisons  in  Fig.  24  arc  in  somewhat  poorer  agreement  than  those  obtained  with  the  simple  (radial) 
Composite  Trend  in  Fig  22a.  However,  the  modified  mixing  mode!  resulted  in  overall  adequate  pre¬ 
dictions  of  the  hydrogen  concentrations.  Clearly,  application  of  results  obtained  with  the  Inverse  Solution 
Technique  could  result  in  a  superior  mixing  model  for  this  type  of  flow. 

The  success  of  this  extremelty  simple  approach,  using  a  constant  value  of  C,  and  less  than 
an  optimum  model)  is  very  surprising.  Obviously,  no  conclusion  concerning  its  general  applicability  can 
be  made  on  the  basis  of  one  case.  Nevertheless,  it  does  suggest  that  this  approach  should  be  investigated 
further:  additional  Cohen  and  Guile  data.  Ref.  14.  are  currently  being  analy/cd.  In  addition,  other  de¬ 
tailed  mixing  and  reacting  flow  data  are  being  sought,  and  will  he  analyzed,  so  that  the  range  of  How 
conditions  and  geometries  considered  may  be  extended.  The  tact  that  the  technique  is  semiem pineal,  so 
that  the  functional  form  of  C  may  be  arbitrarily  selected,  and  made,  for  example,  a  function  of  the  How 
conditions,  suggests  that  it  may  be  sufficiently  flexible  for  successful  application,  to  at  least  limited  ranges 
of  practical  turbulent  reacting  flows.  Of  course,  the  functional  form  of  C  will  vary  greatly  depending  on 
the  form  of  the  mixing  model  selected.  A  model  for  which  C  *  I .  i.e.,  is  independent  of  the  flow  con¬ 
ditions.  would  he  the  ultimate  since  such  a  model  would  adequately  characterize  the  mixing  and  reacting 
How  without  the  need  for  empirical  adjustments. 
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"  Predicted  by  Model  with  *'C"  Factor  of  0.25 
O  —7.0  in,  I 

O  -  10  0  in  I  £*pw,mentil  Values  Calculated  from 

I  Moan  Velocity,  Density  and  Concentration  Diu 
A  —14.0  in.  / 

Composite  Trend 


R  f,Uf7  :3J  ^rnUl  MjSS  Trans,'er  Coefficient.  (.  tnd  Profiles  Predicted  by  the  Model 

1  ’  ‘  :Z!/rT  0i  ?r25  Transition  Re^lon  of  J  Reacting  Hydrogen/Air  Jet. 
(Data  ot  Cohen  and  (,uile.  Rel  14.  High  1  emperature  Vitiated  Air  Case! 
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MODELING  nil  COR  I-  REGION  IN  TURBULENT  AXISYMMI  IRK  QUIESCENT  JUS 
ANDCOI  LOWING  STREAMS 


1  he  core  region  ol  lurhulenC  jets  (  Fir.  25 » is  ol  eons, den. hie  interest  lor  various  practical  aPrl, ca¬ 
tions  For  example,  it  is  in  the  core  legion  where  ignition  ol  reactive  gases  generally  occurs  and  where 
the  greatest  contribution  to  aerodynamic  jet  noise  production  is  made.  In  addition,  the  prediction  ol 
l  ow  held  development  downstream  ol  the  core  region  requires  initial  starting  profiles;  an  estimate  of 
these  profiles  can  usually  he  best  approximated  at  the  point  of  injection,  i.e..  at  the  jet  nozzle  where 
knowledge  ol  the  upstream  conditions  can  lie  used  to  approximate  the  boundary  layer  profiles.  Upon 
specifying  the  initial  protiles.  predictions  of  useful  design  parameters,  o.g..  combustion  efficiency,  can  be 
made  by  solving  the  governing  shear  layer  equations  using  the  methods  such  as  those  reported  in  Ref.  20. 

A  core  model  which  used  a  simple  step-type  slug  profile  for  the  initial  condition  at  the  injection 
station  would  be  useful  because  mean  and  turbulence  parameters  are  rarelv  available  at  the  point  of  m- 
jechon.  However,  turbulence  plays  an  important  role  in  correctly  modeling  turbulent  reacting  flows..e;g 
Kef  lb,  1 7  and  mit-al  conditions  may  significantly  effect  the  turbulence  parameters  of  interest,  eg., 
turbulence  intensity  and  concentration  fluctuations.  Hence,  it  is  possible  that  for  some  flow  conditions 
a  simple  step-type  or  slug  profile  may  be  adequate;  whereas,  for  other  conditions  it  is  necessary  to  consider 
the  details  ol  the  boundary  layer  profiles  at  the  jet  rozzle.  Therefore,  a  twofold  approach  towards  develop¬ 
ing  a  general  core  model  was  taken  in  this  investigation :  results  have  been  reported  in  Refs.  7  and  2 1 

First,  a  model  was  developed.  Ref.  7.  which  started  all  calculations  from  simple  step-type  tslugt 
prolife  and  was  shown  to  give  adequate  predictions  of  mean  values.  A  second  model  was  developed  Ref 
- 1 .  which  used  measured  initial  profiles  and  predicted  both  mean  and  turbulence  parameters  in  the  core 
and  regions  downstream.  A  summary  of  these  models  is  presented  below. 

I  A  CORE  MODEL  USING  INITIAL  STEP-TYPE  PROFILES 


As  shown  in  Ref.  7.  the  generalized  eddy  viscosity  model  was  reasonably  successful  in  the 
ra,l's't'on  rcgl°n  for  w,m’h  it  was  developed.  The  simplest  possible  approach  was  to  assume  a  core  model: 
with  the  same  functional  lorm;  therefore,  as  a  first  approximation,  the  transition  model  was  multiplied 
by  a  constant  factor  (less  than  unity  1  to  correct  for  the  overmixmg  predicted  using  it  in  the  core.  The 
constant  0.4  proved  to  be  satisfactory,  and  the  “core  model"  which  resulted  i  e  e  =  0  4  e 
was  far  more  successful  than  anticipated.  Of  course,  when  applying  for  the  mass  defeat.  Eq.  < 9a? "were' ™ 
changed  Irom  0  to  «(the  extent  of  the  mixing  zone  in  the  transition  region,  to  r,  and  r,  (the extent  of 
the  mixing  /.one  in  the  core).  The  velocity  half  width  value  of  the  eddy  viscosity  at  a  given  axial  station 
in  the  transition  or  similarity  regions  was. 


_  0.018  ( l+2ft(l+e  4  M’>  L’^  /'"  ,,cl  rdr 

ru  +  (D-n  expu  (-0.1 15z) 

Variation  ol  the  eddy  viscosity  in  the  transverse  (radiali  direction  was  given  by 

_L  -  /°u'}  _  1.05-0. 1 5e  '4'6r  ru 

eu  ipu'2  >u  1.0+CL05  <r/ru)7  = 


<9h) 


51 


Bell  Aerospace  Company 


Figure  25.  Schematic  of  Coaxial  Turbulent  Jet  and  Definition  of  Mixing  Regions 
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\  quafion  defined  the  generalized  eddy  viscwsdv  mode  I  developed  lor  the  transitior  region 
I  Ills  model,  together  with  the  assumption  th.il  Sv  j  and  IYj  were  constant,  was  used  tor  predicting  the 
mixing  ot  the  tree  turbulent  shear  l lows  selected  by  the  NASA  l  anglev  Shear  I  ayer  (  onfercnce  Data 
Selection  (  omimltee.  Kel.  7  Success ol  this  model  was  demonstrated  h\  the  reasonable  predictions 
attained  when  iisinp  cither  experimental  or  ship  profiles  at  the  iniedion  Nation.  This  success  suggested 
that  the  genet al  functional  lorut  «>l  the  transition  model  might  be  appropriate  lor  the  core  as  well.  Iyprc.il 
predictions  are  shown  below  tor  a  high  speed  quiescent  let.  ol  Maestrello  and  McDaid  Ret.  22. and  a  ur* 

I  low  mg  hvdrogcn  jet  exhausting  into  a  slower  moving  air  stream.  C  ase  I  A  (  hriss  Ret.  6. 

I  or  the  quiescent  jet.  the  predicted  centerline  velocities  were  at  most  I4r<  greater  than  their 
experimental  values  as  shown  in  Fig.  26.  This  good  agreement  indicates  that  the  initial  slug  velocity  pro¬ 
file  was  adequate  lor  this  case.  However,  as  expected,  even  better  agreement  was  obtained  when  the  actual 
experimental  core  profile  '*'as  used.  These  results  tend  to  validate  the  simple  core  model  for  use  with 
realistic  profiles,  as  well  as  lor  slug  profiles. 

Results  tor  the  high  speed,  subsonic,  collowing  streams  of  hydrogen/air  case  are  shown  in 
Tigs.  27  and  2X.  The  centerline  predictions  obtained  using  both  an  initial  slug  profile  and  experimental, 
core  profiles  underestimated  the  mixing  by  29' l ,  with  significant  difference  observed  between  these  cases. 

C  ah  illations  started  from  the  transition  region  show  slightly  poorer  agreement  was  obtained  4 within  45f'  i. 
The  first  experimental  transition  velocity  profile,  obtained  at  x  D  =  5.34.  is  presented  in  Fig.  28.  C.ood 
agreement  at  this  station  again  demonstrated  the  adequacy  ot  the  simple  core  model.  However,  these 
results  show  that  the  eddy  viscosity  model  does  not  include  all  the  complexities  required  for  predictions 
of  all  core  f  low  s. 


2  PI  VISIT)  H)!)Y  VISCOSITY  AND  TURBULI  NC  F  INTENSITY  MODFLS  FOR  Till  C  ORF 
,  IT, ION 


As  part  o|  a  PhD  thesis.  Ret.  21.  the  eddy  viscosity  models  of  Refs,  7  and  13  were  modified 
in  an  attempt  to  improve  predictions  for  the  transition  region  of  quiescent  jets,  wakes,  and  cotTowing 
streams:  this  model  also  was  generalized  to  model  the  core  region.  The  modified  model  contained  few^er 
empirical  functions  than  the  previous  model,  i.e..  it  did  not  contain  an  empirical  expression  for  the  char- 
t  or  i  si  ic  length,  or  a  I  unction  to  account  tor  flows  in  which  the  velocity  of  the  external  stream  is  greater 
than  the  jet  velocity.  1  hese  modifications  have  led  to  an  improved  ability  to  make  predictions  of  tur¬ 
bulence  intensitx  in  the  core  region  for  several  cases.  Application  of  this  model  to  the  core  is  discussed 
below. 


The  core  may  he  divided  into  three  distinct  regions  as  shown  in  Fig.  29.  The  initial  reeir^u- 
lation  /one  i  Region  1 1  tor  most  jet  mixing  problems  is  smalt  and  is  assumed  to  have  negligible  effect  on 
Mie  downstie.ini  flow  development  Flows  in  which  recirculation  is  important  require  analysis  usinsi  r he 
Nav ter-Sirokes  equations  since  both  longitudinal  and  transverse  pressure  gradients  must  be  considered 
Such  an  anaKsis  was  box  and  the  scope  of  this  investigation  but  is  recommended  as  a  subject  fora  future 
investigation.  The  problem  of  modeling  the  core  region  thus  reduces  to  modeling  the  quasi-two-dimeriMonal 
shear  layer  ( Region  II  i  and  the  region  downstream  (  Region  HI  i  where  jxbymnietry  effects  become 
im  porta  nt. 
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Figure  28.  Predicted  and  Kxperimental  Axial  Velocity  versus  Distance  from  Centerline  at  7.1  D  =  5.34. 

Data  of  Chriss,Ret.  6.  Case  1 A 
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Figure  29.  Core  Region  of  Axisymmetric  Turbulent  Jet 
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I  lie  incompressible  c  descent  jo?  data  reported  by  Sami,  Rot.  23.  was  chosen  as  tho  starting 
point  tor  tins  core  model  development  since  it  is  tho  most  detailed  core  data  available.  I  he  result  mo  edd\ 
Mscositv  amt  turbulence  intensity  models  below.  Iqs  10  and  I  I.  were  then  applied  to  compressible  ipues 
vent  and  coflowing  jet  data  and  touud  to  give  reasonably  jv  curate  predictions  ol  the  mean  vclocitv  Ik  Id 
iti  addition  Jo  accurately  correlating  the  shear  stress  and  turbulence  intensity  reported  by  Sami,  as  shown 
in  F:ig  30  and  3 1 . 
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A  comparison  of  experimental  data  and  the  predictions  tor  a  high-speed  quiescent  jet.  Maestrello 
and  McDaid.  Ref  22.  in  l:ig.  32.  Agreement  is  quite  good:  predictions  for  this  case  are  within  16  percent 
"t  the  data  (  omparisnn  t oi  a  cottoning  h\dr<>gen  jet  exhausting  into  a  slower  moving  air  stream.  (.  nriss 
(  ,:nc  i  \  tiie  data  used  1<  >r  use  Client  nf  the  step-type  initial  profiles.  is  shown  in  lies.  33  and  34.  I  hese 
ligures  show  that  predictions  lor  this  cot  lowing  hydrogen  air  data  of  (  hriss  are  witlun  33  percent  of  the 
experimental  centerline  velocities,  which  appears  adequate  for  main  hardware  applications. 
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O  Exp,  Sami 
11  Predicted 


Figure  30.  Predicted  and  Experimental  Turbulent  Shear  Stress  at  Velocity  Half  Width 
for  Data  of  Sami,  Reference  23 
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Figure  3 1.  Predicted  and  Pxperimental  Axial  Turbulence  Intensity  Profiles 
at  2  =  3  and  10  for  Data  of  Sami,  Reference  23 
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Figure  34.  Mean  Velocity  Profiles  at  r  D  =  6.7  for  Data  of  Chriss.  Ref.  6,  Case  t  A 
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(  TWO-PHASI  MIXING  STl'DY 

T,w  benefits  tor  rocket  and  ramjet  engine  performance  that  nuv  he  derived  Iron,  the  injection  of  . 
slream  ol  -did  particles  into  a  secondary  stream  can  he  assessed  only  alter  the  complicated  mixing  and 
combus  ton  processes  are  reasonably  well  understood.  This  part  of  the  program  deals  specifically  with 
the  tmiu  ilv  numicjl  aspects  ot  the  mixing  process. 

Studies  elsewhere  apparently  have  Iven  limited  to  particle  injections  in  the  direction  of  the  average 
low  or  at  right  angles  to  it  <  Ret.  -'41  In  particular,  the  injection  region  has  not  been  considered  in  detail 
In  a  practical  combustor,  injection  of  fuel  particles  generally  would  be  at  some  intermediate  angle  to  the 
primary  stream  to  maximize  mixing  and  minimize  losses.  Initially,  the  particles  tend  to  penetrate  into 
the  primary  t  ow  as  a  result  of  their  inertia.  As  the  particle  trajectories  become  more  or  less  parallel  to  the 
;!,f,lV''°n  ,ur,hcr  n"''nf  's  controlled  mainly  by  the  slower  process  of  turbulent  diffusion: 
this  tar  tick!  behav tor  ol  the  particles  is  currently  being  studied  elsewhere  ( Ref.  25 1.  I  he  investigation 
discussed  here  is  concerned  mainly  with  the  region  near  the  injection  point,  where  the  particle  acceleration 

it!1' mn  rV  V'VT  *(ranl‘adS  ,0  a  markcd  curva,ure  ,ho  trajectory.  This  injection  region  appears  to  have 
Hk  most  important  influence  on  overall  combustor  performance  (Ref.  8). 

T1,°  c;ornplete  problem  involving  particle  slip  with  respect  to  the  gas.  ignition  delay,  and  combustion 
as  sll  as  turc  ilcnt  mixing  ol  the  gases  and  ol  the  particles,  is  extremely  difficult  to  solve.  Therefore  il 

oneT.Te  5f  ‘  T  r"h ^7  f'°W  ",0dels  which  consider  onl>  elected  aspects  of  the  problem  at 
one  time.  The  results  ol  such  studies  should  provide  an  insight  into  the  relative  importance  of  the  processes 


genti  ilh  nt,  a  fl  .  ,  Pm  ,  Pr7rUm  dcals  wi,h  ,he  dVnan,i^  ol  Par,|cles  injected  nontan- 

gentially  into  a  flow  Several  problems  related  to  the  behavior  of  single  particles  were  investigated  The 

major  aspects  ol  this  work  was  described  in  the  proposal  for  the  current  contract.  Ref.  26;  a  final  report 

l  “  ,,ci"g  p:cpared  AI,hou*h  the  of  Single  particles  gives  valuable  insigl.t,  it  is 

important  to  investigate  also  the  injection  of  a  stream  of  particles. 

tns„  lcTaPr  l’’  two-dimensional  gas-particle  jet  was  represented  by  a  pure  gas  jet  of  appropriate 

uKomnr  .ss'ih  CU  '  "  ''7?  C°n,’"ed  para"cl  Wa,,s  ,Fi*  3?a»  The  "ow  "as  treated  as 

whe«K  .h  •  nivt  r  a  K‘  pr°b,eni  was  at,acked  bv  the  OH-thod  of  conformal  transformations. 

I  n  , ,  f.Pj  P  "  'S  ,,anS'°rnlcd  1,1,0  ;1  P'ane  in  wh,dl  'he  coordinates  represent  equipotential 
at  s  rum  in  n  am  stream  and  the  jet  then  torni  two  separate  flow  fields  (Fig  A Sb t  and  the 

proper  boundary  ,  uions  mist  be  satisfied  along  the  walls  and  jet  boundaries.  A  computer  pro-ram 

wav  prepared  and  al  e  results  have  t'een  compared  with  experimental  observations  I  see  Appendix  (  )  i 
report  on  this  etfort  w  •!■  be  issued.  Ret.  2".  1 

To  provide  adequate  residence  time  for  particles  to  ignite  and  burn  in  a  combustion  chamber  it 
appears  promising  to  inject  the  particles  with  a  velocity  component  opposing  the  main  stream  m  the  com- 
us  ion  e  umber,  tut  one  must  anticipate  that  such  a  scheme  would  imply  certain  penalties  The  problem 
requires  the  analysis  ol  simultaneous  mass,  momentum,  and  heat  addition  to  a  flow  and  represents  an 
extension  ol  earlier  work  on  pure  heal  addition  (Refs.  28  and  29) 
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Figure  35a.  Two-Dimensional  Jet  Injection 


Mam  Flow 


Jet  Flow 


Figure  35b.  Transformed  Plan  of  Figure  35a 
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Three  different  assumptions  lor  the  main  How  were  considered:  (a)  constant  subsonic  How  and 
reservoir  temperature.  tbi  supersonic  How;  and  (cl  constant  reservoir  pressure  and  temperature.  In  all 
cases,  the  flow  from  the  combustion  chamber  was  assumed  to  ho  completely  mixed  and  burned  and  to 
exhaust  directly  into  the  atmosphere.  The  flow  variables  were  determined  for  various  conditions  ot  the 
injected  fuel  jet.  The  three  cases  are  indicated  in  Figure  3<>.  where  M0  and  in0  represent  the  Mach  number 
and  the  main  How  rate  without  fuel  injection,  while  Tr  and  pr  represent  the  reservoir  temperature  and 
pressure  of  the  main  flow.  After  injection  of  the  jet  How  nvj  =  Jiii0,  the  Mach  numbers/ upstream  and 
downstream  of  'he  combustion  /one  are  denoted  by  M,  and  M: .  The  velocity  of  the  fuel  jet  is  given  by 
Ujar  where  af  is  defined  as  the  reservoir  speed  of  the  sound  of  the  main  How.  The  values  of  Uj  may  Ik* 
positive,  zero  or  negative  depending  on  whether  injection  is  in  the  downstream,  lateral  or  upstream  direc¬ 
tion.  All  three  cases  in  Figure  36  are  shown  for  negative  values  of  Uj.  The  fuel  jet  has  a  reservoir  tempera¬ 
ture  aTr  and  the  amount  of  heat  released  by  combustion  per  unit  mass  of  injected  jet  is  denoted  by  q 
which  can  be  expressed  in  terms  of  the  enthalpy  of  the  main  How  by  h  =  q/c^T^  The  pressure  in  the 
external  discharge  region  is  pe>  and  the  boundary  conditions  that  must  be  satisfied  at  the  discharge  exit 
of  the  duct  are  given  by 

M2  <  1  andp:  =  pc  or  M2  =1  and  pc  >  pe 

Conserv  ation  of  mass,  momentum,  and  energy  then  provides  a  sufficient  number  of  equations  to  solve  for 
all  How  conditions. 

Case  fa)  represnts  the  condition  in  which  the  main  flow  to  the  combustion  chamber  is  supplied 
from  a  supersonic  diffuser  te.g..  supersonic  ramjet).  The  analysis  then  indicates  that  the  thrust  is  indepen¬ 
dent  of  the  direction  of  the  injected  fuel  jet.  Thus,  the  negative  thrust  contribution  of  an  upstream-directed 
fuel  jet  is  exactly  compensated  by  a  reduced  shock  loss  in  the  diffuser. 

Case  (b)  represents  injection  into  a  supersonic  flow  (e  g.,  scranijet).  As  long  as  the  combined  effect 
of  mass,  momentum,  and  heat  addition  does  not  exceed  a  critical  value,  the  upstream  flow  remains  con¬ 
stant,  and  the  thrust  of  the  model  considered  here  is  given  by  the  sum  of  the  thrust  without  injection  and 
that  of  the  injected  jet  regardless  of  the  amount  of  heat  released.  Critical  injection,  at  which  an  upstream- 
traveling  shock  wave  is  formed  to  make  the  upstream  velocity  subsonic,  is  determined  by  a  relationship 
between  the  parameters  J.  Uj  and  a  +  h.  Figure  37  shows  the  critical  values  Jcrjt  for  lateral  jet  injection 
( U;  =  0)  as  function  of  the  upstream  Mach  number  for  several  values  of  a  +  h  and  for  a  gas  with  a  specific- 
heat  ratio  7  =  1.4.  If  Uj  is  not  equal  to  zero,  these  values  of  J crj j  are  further  modified  to  Jcrjt  +  as 
indicated  by  Fig.  38. 

Case  id  reprevenK  injection  into  a  flow  from  a  constant  pressure  reservoir  (e  g.,  a  stationary  burner 
in  a  power  dationi  with  pr  p  ~  K.  In  this  cuse.  the  thrust  of  the  system  depend^  in  a  complicated  manner 
on  the  five  parameters  K.  J,  Uj,  a  +  h,  and  y.  A  computer  program  was  prepared  to  explore  the  influence 
oi  ! he m.’  quantities,  and  some  of  the  results  are  shown  in  Figs.  39.  40,  and  41.  Let  v  denote  the  ratio  of 
the  thrust  for  some  value  of  Uj  to  the  thrust  for  lateral  injection  (Uj  -  0).  Figure  39  shows  ^  as  a  function 
of  the  compression  ratio  K  for  several  values  of  J  and  Uj  and  for  h  =  1 0.  a  =  1 .  and  y  =  1.4.  The  effect  of 
jet  direction  is  indicated  in  Fig.  40  and  is  almost,  but  not  exactly,  symmetrical  with  respect  to  the  sign  of 
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Figure  36.  Mass,  Momentum  and  Heat  Addition  to  a  Gas  Flow 
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Figure  39.  Fitted  of  Compression  Ratio  on  Thrust  Ratio  $  =  Thrust/Thrust  for  U: 
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Figure  40.  Effect  of  Jet  Direction  on  Thrust 
0t  =  10, a  =  l.r=  1.4,  J  =  0.1 ) 
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Figure  4 1 .  Separate  Effects  of  Heat,  Mass  and  Momentum  Addition  (  y  =  1 .4. 
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t  is  interesting  to  compare  the  mass  How,  ax  affected  by  mass,  momentum  and  heat  addition 
separately,  Figure  4  shows  the  mass  How  ratio  m,  m„  as  a  function  of  the  compression  rat.o  K  for 
„  .  '  '"H  J _4n  r,K;  "'Z  represent  ...ass  add.hon  U-  O  il  for  the  three  orientations  of  the 

r  ,"  Jo°r  .  Ihe  broken  line  at  the  center  represents  the  effect  of  heat  addition  alone 
lit  *  ,r0C  n0"'  CUrVCS  ,ndica,c  ,hc‘  co"'bl"«» «?»»«»  »l  mass,  momentum,  and  heat 


The  three  basic  cases  considered  are  those  of  mam  importance  for  combustor  design.  A  detailed 
descr.pt, on  of  .he  method  of  analysis  will  be  given  in  a  paper  which  is  being  prepared  Ref  30. 

Flow  and  combustion  instabilities  in  rockets  and  other  devices  in  which  gas-particle  mixtures  are 
employed  often  are  sustained  by  the  feedback  provided  bv  pressure  waves  in  the  system,  [t  is  therefore 

3S1 VTTtU  ,hC  PrOPCr,iCS  ir*  WaVCS-  am*  a  SUrVCy  Pa'"r  thiS  ™  -  P-P-d  ^r  • 

publication  Kef.  3 !  This  survey  covers  all  types  of  pres»  »  re  waves: 

of  hieher'nr  ‘’‘"Si*  Wi''Ch  produ'*  a  transi,ion  ft™"  one  now  condition  to  another 

b[  a„  X,  „  .  ,CSC  3  Prat,ica,,y  disc '"•tin  nous  pressure  rise  at  the  wave  front  followed 

,K-r*d  shoektr  3  l0n  Z°ne  °r  Pr°dl,CC  the  en,ire  ‘  'i,ion  wi,ho"‘  discontinuity  in  the  form  of  d,s- 

(2)  Sound  waves  which  represent  periodic  disturbances  of  small  amplitude  and  fixed  frequency: 
pressure'  3'  L3rge'3n,pNtU<le  Waves  of  arbitrary  wave  form  which  may  produce  a  net  rise  or  drop  of  the 

The  analysis  of  fuel  jets  injected  in  the  upstream  direction  discussed  in  the  preceding  section  does 

intecUd  solid  particles.  Large  particles  penetrate  further  upstream  but  also  take  longer  to  burn  than  sm  ,11 

o  ("t  he  co  n !  bust  ion  di  am  be  !°  7n  l  h*- '  TtT™  PM'  ""gh'  exist th:,t  leads  t0  a  ™n™um  length 

h ...  .  nt?  m  u  ^  P3rt  °‘  'he  Pr°gram-  deta,ls  ot'  combustion  will  not  be  considered 

but  a  plausible  allowance  for  the  burning  time  may  be  made  based  on  available  data  Refs.  32  to  3S. 

face  of  ,hi  ™ir.T  01  ParI"dCS  'S  ,limitcd  by  ,he  ra,°  Jl  wbich  ^seous  reactants  can  diffuse  to,  the  sur- 
fact  ol  he  particle.  In  pnnetple.  ,t  ,s  therefore  possible  lo  increase  the  burning  rate  by  letting  the  particles 
sl-P  with  respect  .0  the  surrounding.  In  tins  manner  the  surface  o  .ontmuouslv  exposed  to  a  hmher  con 

to"  tlm’?  °,i  T  rearan!  ?3S'  Und  ,lK'  d,,U,sion  limi,a,io-1  ls  P‘"»>  overcome.  An  attempt  is  bring  made 
!  n  m  m  r  r"  vlCiT  ,erived  lr°m  par,,de  ^  and  <"  <M’l<-re  techniques  tor  produemu  such 

:l:i;::;::;:fo;::,;;;^uu;y  zt . .  ,s  < . . — - -  *  ^  ,1,. 

n  MODELING  Of  Mb  AN  VALL'F  DATA 

I  CENTERLINE  DECAY  EXPONENTS 

Universal  jet  wulth  growth  and  centerline  decas  laws  are  known  lor  constant  density  smelt* 
stream  jets  and  wakes  the  similarity  region.  Reference  36.  It  would  be  useful  if  comparable  exprewons 
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we 10  available  tor  iihmo  complex  flows  mkIi  .is  variable  densitv  two  stream  tets  I  horotoro  methods  proposed 
to  model  centerline  mass  traction  decav  in  turbulent  axisymmctric  lets  wore  evaluated  and  arc  discussed 
bo  low 

Zakkay.et  al  Kef  37  have  suggested  that  the  centerline  mess  traction.  Y^.  may  be  determined 
front  the  relation 

Y  .  *  (/,/  Xm  . 

C  C 

wln-rc  in  =  2.0  and  /  iv  the  potential  core  length  and  /  is  avial  distance  downstream  of  the  injection  station. 
Abramovich  et  al  Ref.  3X  have  suggested  the  same  functional  form  as  I  q.  ( 1 2 1  but  state  nt  =  I J  is  the  correct 
value.  Schet/  Ref.  3'»  has  suggested  that  the  decav  exponent  can  be  correlated  with  the  ratio  of  the  injected 
mass  flux. 

( pi’ )}  i  pi1  te  where  he  showed 

m  =  2.0  loriplMj  fplh  <1.0 
and 

m  =  1.0  foripTij  (pUie  >  10 

Cohen  and  Guile  Ref  I  d  state  that  the  decav  exponent  lor  mass  fraction  is  unity  and  independent  of  the 
den m! \  ratio,  p  pc  if  the  \elocity  ratio  l  j  l’e  is  near  unity.  Their  results  suggest  m  =  1 .0  for  0.67  m  L’j-L’0 
<2.Ch  ' 

Table  4  lists  the  values  of  ipUk  lpU>c.  Uj  t:c  and  in  for  the  data  used  by  Schet/.  Ret  MK  and 
Cohen  and  Guile.  Ref.  14.  in  arriving  at  tneir  conclusions.  For  completeness,  data  not  considered  by  Schet/ 
or  Cohen  and  Guile  also  was  analyzed,  i.e..  the  data  of  Forstall.  Ref.  40.  Zawacki  and  Weinstein.  Ref  l  I. 
Ragsdale  et  al.  Re*.  41.  IVSou/a  et  al.  Ref.  42.  Chriss.  Ref.  6.  Fggers  and  Torrence,  Ref.  5.  and  Keagy  and 
Weller.  Ref.  43.  \  total  of  64  different  cases  are  summarized  in  Table  4  from  nine  different  investigations 

The  centerline  mass  fraction  decay  exponent  was  obtained  by  using  a  least  squares  fit  ot  I  q  ( 12) 
and  vp!  mg  for  /  and  m .  see  Ref  21  for  details.  Thi-  technique  differed  from  past  methods  of  obtaining  /  . 
.snd  Di  where  a  straight  one  w  p  mi1  p}\  drawn  through  data  on  a  iog-|*»e  plot,  and  its  slope  deter n  med  graph* 
ie.:*i\.  t  one  the  k  »si  square"  tu  technique  is  important  since  the  difference  between  a  decay  exponent  id  m- 
1  5  and  2.0  is  only  an  angle  of  ",c  when  plotted  on  log-log  paper.  A  judgment  as  to  when  to  neglect  points 
near  the  core  region  must  he  made  when  making  these  plots.  I  his  judgment  is  made  in  a  more  quantitative 
manner  using  the  least  squares  lit  technique  since  it  is  required  that  the  calculated  values  ot  zL.  and  m  give  a 
residua!  of  less  than  0.01 ;  therefore,  points  near  the  core  were  automatically  deleted  when  including  them 
resulted  in  residuals  greater  than  0.0 1. 

F xarninat ion  of  the  \ allies  of  m  showed  that  these  values  change  significantly  depending  on  the 
type  of  flow  considered.  Clear !>  universal  values  of  m  =  2.0.  Zakkav  et  al.  Ref.  37.  or  m  =  1 .7  Abramov T 1; 
et  al.  Ref.  38.  are  not  supported  by  these  results.  The  value  of  m  appears  to  be  dependent  on  the  ratios  fpl  >j 

\\%  and  m  other  parameters  as  well.  In  addition,  the  eti-ect  of  the  expermient.il  .on  turn*' 

ration  which  also  n  a>  pia\  a  role  m  determining  the  value  ot  m.  w as  discussed  by  Abramov ich  el  al.  Ret  e* 
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Figure  4 2  shows  the  values  o  n  plotted  versus  the  mass  f  lux  ratio  < pi’ >(  (plMc  aiul  the  fiend 
observed  by  Schet/.  Ret  W  I  or  clarity.  jet  like  flows  iU  I ^  >  1 .01  are  represented  by  open  symbols 
whereas  wjkehke  flows  d  (  l  ^  •  -  I 1  are  denoted  as  solid  symbols  It  only  the  jet  like  flows  are  considered, 
this  plot  shows  that  the  trend  observed  by  Sehet/  is  reasonable  except  that  a  more  representative  value  oi 
in  tor  (pi  )}  ipl’)0  <  I  appears  to  be  nearer  1 .7  Ilian  2  0.  as  suggested  by  Abramovich  et  at,  Ref  3N  On  the 
other  hand  it  both  jet  and  wake  like  data  are  considered,  no  discernible  correlation  between  in  and  (pi  »  f 

Ih.*  obtained  Interestingly.  considering  only  wake*ltke  data  (solid  symbols)  suggests  (hat  maximum 
values  of  m  may  occur  near  (pl’f  (pi  *  1.0.  However,  the  scatter  tti  Fig.  42  precludes  any  definite  con¬ 
clusions  and  strongly  supports  the  recommendation  that  additional  experimental  data  be  obtained,  using 
advance  diagnostic  techniques,  such  as  the  laser-Doppler  system. 

Figure  43  shows  the  values  of  m  plotted  versus  the  velocity  ratio  VyVc  and  the  trend  observed 
by  Cohen  and  Guile.  Ref.  14  (  ases  where  p-pe  >  1.0  are  represented  by  open  symbols,  whereas  cases  with 
Pj  Pc  ^  I  0  are  designated  by  dosed  symbols.  Note,  these  results  do  not  support  the  observation  ot  (  often 
aiid  Guile  that  m  =  !  over  the  range  O.o7  ^  ^  2.0.  In  trying  to  explain  this  discrepancy,  a  detailed 

examination  ot  the  decay  exponents  by  (  ohen  and  (mile  was  made  and  compared  with  values  obtained  in 
this  study. 


C  ohen  and  (mile.  Ret.  14.  have  obtained  reacting  hydrogen  air  turbulent  jet  data  for  four  dif¬ 
ferent  sets  of  tree-stream  conditions.  The  tree-stream  total  temperature  was  approximately  1750°K  in  two 
^ases  and  I  5b0  K  in  the  other  two.  The  density  and  velocity  ratios  of  each  cast,*  were  kept  approximately 
constant.  Figure  44  shows  the  centerline  mass  fraction  decay  for  these  four  different  cases.  In  Fig.  44.  “high 
NV'  refers  to  a  case  with  the  higher  total  temperature  ( 1 750'Kt  where  the  vitiated  external  stream's  ()2  de¬ 
pletion  was  lepljeed  by  \2.  “High  air  signifies  that  the  depleted  02  was  replenished  with  make-up  oxygen 
and  the  total  temperature  also  was  1 750°K.  Note.  Fig.  44  shows  that  the  Cohen  and  (mile  data,  as  well  as 
the  data  of  Alpimeri.  Ref.  41.  is  either  in  the  core  or  just  downstream  of  the  core  region,  i.e..  in  the  transition 
region.  However,  meaningfully  asvmpototic  decay  exponents  cannot  be  defined  in  the  transition  region  where 
m  ranges  from  zero,  just  at  the  end  ot  the  core,  to  some  finite  value  further  downstream,  therefore,  the  decay 
exponents  obtained  by  Cohen  and  Guile  are  not  asvmpototic  decay  exponents. 


The  uncertainty  involved  in  calculating  the  decay  exponent  of  the  centerline  mass  fraction  as 
demonstrated  from  the  above  discussion  leads  to  the  following  conclusion.  The  decay  exponent  of  the  center¬ 
line  mass  fraction  cannot  be  correlated  with  the  velocity  ratio  V  1/  or  the  mass  flux  ratio  ipUl  /  (pc  j 
aione  using  existing  experimental  data.  I(  such  a  correlation  exists,  data  from  the  far  downstream  region 
similarity  region!  covering  a  wide  range  of  velocity  and  density  ratios  must  be  used.  Two  facts  which  make 
o;c:;  data  d;M icult  to  oht4,m  using  conventional  sampling  techniques  are-  I  1 1  in  the  far  downstream  regions 
the  par.imvk  t<  being  measured  arc  small  and  hence  subject  to  large  experimental  uncertainties,  and  (  2i  the 
simulation  ot  an  infinite  external  stream  becomes  increasingly  difficult  far  downstream  due  to  external  in¬ 
fluences.  e  g.,  walk  Therefore,  to  obtain  such  data,  superior  diagnostic  techniques  like  hot-wire  and  laser- 
Doppler  anemometers,  will  he  required  as  well  as  large  facilities. 


2  RFICHARDTS  INDIX TIVF  THEORY 


Reichardt  (see  Hin/e.  Ref.  36>  noted  that  velocity,  temperature  and  concentration  distributions 
across  mixing  /ones  closely  follow  Gaussian  error  functions.  By  assuming  that  a  turbulent  transport  pr-vess 
is  a  statistical  process  exactly  analogous  to  molecular  transport,  it  follows  that  the  differential  equation  for. 
say.  momentum  rtuM  be  identical  to  the  heat  conduction  equation.  Ilin/e  Ref  3b  has  shown  that  Reichardt's 
liypnfl ,,-sis  conflicts  with  Newton''  relativity,  which  states  that  torses  m  a  meJumca!  system  must  tv  m- 
ds  ps  ndent  ot  the  addition  ot  constant  velocity.  Nevertheless.  Reichardt  s  theory  is  appealing  since  is  Jesuits 
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Figure  43.  Centerline  Mass  Fraction  Decay  F.xponent  versus 
Velocity  Ratio  Uj'Ue 
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m  line.  *  equations  and  din's  not  require  similarity  assumptions  nor  alterations  lo  account  lor  compressibility. 
Xrul  niovt  importantlx.  it  had  given  eorreet  agreement  with  data  lor  single -tree  jets  mixing  with  a  quiescent 
ambient 


1  rent aeoste  and  Slor/j.  Ret.  45.  have  suggested  I -wilier  pursuing  the  Reiehardt  theory  to  pro 
hiems  involving  mass,  momentum,  and  energy  transport.  In  the  study  reported  herein,  two  types  of  initial 
How  conditions  were  found,  for  which  it  appears  that  Reiehardt's  theory  predicts  trends  inconsistent  with 
the  reported  data.  However,  the  results  shown  he  low  are  somewhat  inconclusive  since  pressure  gradients 
in  the  near  core  region  may  account,  at  least  in  part,  foe  this  contradiction  of  the  Reiehardt  theory.  I  valua¬ 
tion  of  the  importance  ot  pressure  gradients  could  not  be  determined  Irom  the  reported  experimental  data. 

Figure  45  shows  the  dimensionless  momentum  flux  as  obtained  from  the  analysis  of  the  data  of 
Chriss.  Ref.  6.  These  results  show  that  for  the  ease  oi  a  light  gas  Ul2 )  injected  into  a  slower  moving  heavy  gas 
tairl  the  momentum  flux  tpU  ■  *c  does  not  exhibit  a  C.aussian  behavior.  Peters.  Ref.  46,  also  has  noted  a 
trend  differing  from  Gaussian  behavior  lor  this  same  data.  This  result  may  be  explained  by  noting  that  the 
dciiMtv  on  the  centerline  is  increasing;  whereas  the  centerline  velocity  is  decreasing,  resulting  in  a  nomnono- 
to.tu  behavior  of  pi* ,  A  similar  trend  was  observed  for  the  data  of  Zawacki  and  Weinstein.  Ref.  I  I .  in 
which  a  heavy  gas  (lreon-1  2)  was  injected  into  a  faster  moving  light  gas  (air). 

It  should  he  noted  that  Reiehardt'*  theory  considers  the  time-averaged  value  of  the  momentum 
flux  rather  than  the  product  of  mean  density  and  mean  axial  velocity  squared  as  was  the  case  in  Fig.  45. 
Suuepl  *  pi  *  +  p  ir  +  2Lp  u  +  pir .  the  last  three  terms  of  this  equation  were  not  considered  in  the 
above  discussion.  Neglecting  these  terms  appears  to  be  reasonable  for  the  data  considered,  and  is  consistent 
with  the  boundary -laser  approximations. 


80 


VI 


Dimensionless  Momentum  Flux  on  Centerline  lor  HydrogeiCAir  Jet 
Data  of  Chris*.  Reference  6 
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APPENDIX  A 

I  URBULl  NT  MASS  TRANSFER  (OI  IMCII  NTS  OHTAINM) 
USING  INVERSE  SOLUTION  TK  HNIOUF 


Section  III-A-3  discussed  the  technique  tor  the  direct  determination  of  turbulent  transport  coefficients. 
It  also  presented  typical  results.  Here,  Figures  A- 1  through  A-IO  present  the  transport  coefficients  obtained 
lor  all  10  cases  summarized  in  Figure  13.  Also  presented  are  the  integral  mass  and  momentum  balances  ob¬ 
tained  at  each  axial  station  for  which  data  were  available.  These  balances  permit  a  critical  assessment  of  the 
consistency  of  the  experimental  data  so  that  the  degree  of  confidence  that  can  be  associated  with  the  given 
set  of  transport  coefficients,  may  be  evaluated. 

Figures  A- 1  through  A-10*  are  plots  of  the  raw  computer  output  resulting  from  direct  application  of 
the  Inverse  Solution  Technique.  In  each  case,  the  figure  designated  “a**  is  a  plot  of  the  mass  transfer  coef¬ 
ficient.  £.  and  the  figure  designated  “b"  is  a  plot  of  the  eddy  viscosity  that  is.  the  momentum  transfer  co¬ 
efficient  >.  (  A  hea\y  line  was  drawn  in  each  figure  to  represent  the  most  significant  variations  of  the  mass 
transport  coefficient,  while  ignoring  less  significant  trends  and  those  effects,  discussed  in  References  H  and  0. 
demonstrated  to  yield  spurious  results.  The  most  important  of  these  effects  are: 

•  The  Abnormally  high  nr  low  values  of  the  transport  coefficients  are  often  obtained  at  the  initial 
and  final  axial  stations  which  results  from  the  difficulty  in  fitting  and  differentiating  experi¬ 
mental  data  ai  the  end  j  oints.  and 

•  Regions  in  the  vicinity  of  the  centerline  and  the  free  stream  in  which  radial  concentration  and 
velocity  derivatives  approach  zero,  may  result  in  anomalously  high  transport  coefficients  as 
indicated  by  Eq.  (A- 1  >, 

.  _  % 

3U/3r  ’  *  3  Y/dr  <A'!) 

where  t(  and  r^  represents  the  Reynolds  transport  terms  for  the  turbulent  transport  of  momentum 
(Reynolds  stress)  and  mass  respectively;  Reference  1 5. 

•  Symmetry  conditions  require  that  r(  and  r^.  and  radial  first  derivatives  of  Y,  U.  Ts.  and  p  be 
zero  at  the  iet  centerline:  therefore,  values  of  the  transport  coefficients  alone  the  centerline  can¬ 
not  he  obtained  directly  using  I  n.  (A-!  i.  However,  transpou  coefficients  at  the  centerline  may 
be  extrapolated  Jrom  results  obtained  in  its  vicimu  or  from  radi.d  second  den\.m\es  at  the 
centerline  t  Reference  ch  Untortunateh  .  second  derivatives  p. '  .rally  are  questionable  when 
obtained  from  numerical  differentiation  of  experimental  date,  -  >  That  extrapolation  appears 
preferable. 


The  Composite  Trends  plotted  in  Figure  13  were  obtained  from  Figures  A-l  through  A-10.  For  sim¬ 
plicity  and  because  the  influence  of  momentum  transfer  was  shown  in  Section  lll-A-2  to  be  far  less  significant 
than  mass  transfer,  the  same  Composite  Trend  was  used  to  represent  both  £  and  e  that  is.  Scj  was  arbitrar¬ 
ily  assumed  to  be  unity  for  the  Composite  e.  In  some  cases,  especially  those  of  Zawaeki  and  Weinstein.  Fig¬ 
ures  7  and  8.  very  little  judgment  was  required  to  draw  the  Composite  Trend;  in  other  cases  it  was  selected 
rather  arbitrarily.  For  the  very  low  speed  flows  in  which  jet  velocities  were  le»s  than  5  ft  'see.  the  intcer.il 
mass  and  momentum  balances  generally  were  extremely  consistent:  for  example  the  momentum  mtceral  in 


*In  these  Figures  z  -  z/Dj ,  where  Dj  is  given  in  Table  I 
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Ma»  Transfer  Coefficient  Profiles  Determined  by  Differentiation 
xpenmental  Data.  ( iinss.  Ret.  0.  Air  ilydrogen  M.xing  ( I  A» 
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Figure  Alb.  Eddy  Viscosity  Profiles  Determined  hv  Differentiation 
of  Experimental  Data.  C'hriss  Ref.  6,  Air; Hydrogen  Mixing  ( 1  A) 
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f  igure  A2a.  Mass  Transfer  Coefficient  Profiles  Determined  by  Differentiation 
ol  Experimental  Data.  Cliriss.  Ref.  6.  Air/I  lydrogen  Mixing ( IB) 


Figure  A 2b.  Mass  Viscosity  Profiles  Determined  by  Differentiation 
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Figure  A3b.  Eddy  Viscosity  Profiles  Determined  by  Differentiation 

of  Experimental  Data.  Oiriss,  Ref.  6.  Air/ Hydrogen  Mixing  ( lO 
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Figure  A4b.  Eddy  Viscosity  Profiles  Determined  by  Differentiation 
of  Experimental  Data.  Cnriss,  Ref.  6,  Air/Hydr^gen  Mixing  ( ID) 
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Figure  A5o.  Eddy  Viscosity  Profiles  Determined  by  Differentiation  of 
Experimental  Data.  Chriss,  Ref.  6,  Air/Hydrogen  Mixing (II:) 
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figure  A5d.  fddy  Viscosity  Profiles  Determined  by  Differentiation  of 
f xperimental  Data.  Chriss,  Ref.  U  Air/Hydrogen  Mixing  (It) 
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Figure  A6b.  Eddy  Viscosity  Profiles  Determined  by  Differentiation  of  Experimental  Data. 
Zawacki  and  Weinstein,  Ref.  1 1 .  Air  Freon  Mixing  <  204F) 
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f  Viscosity  Profiles  Determined  by  Differentiation  of  Experimental  Data 
Zawaeki  and  Weinstein.  Ref.  I  1,  \ ir  Freon  Mixing  (20SF) 
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Figure  A*..  Mass  Transfer  Coe.fie.ent  Profiles  Determined  by  Differentiation  of  Experimental  Data 
Zawacki  and  Weinstein,  Kef.  1 1,  Air/Freon  Mixing  (206F) 
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Liiirc  ASh.  Eddy  Viscosity  Prolilcs  Determined  by  Differentiation  of  Experimental  Data 
Zawacki  and  Weinstein,  Ref.  1  I,  Air/Freon  Mixing  ( 206F) 
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Impure  A9b.  Eddy  Viscosity  Profiles  Determined  by  Differentiation  of  Experimental  Data 
Eggers  and  Torrence,  Ref.  5,  Air  Air  Ethane  Mixing  (Case  1  ) 
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Figure  7  was  identical  at  the  first  three  axial  stations,  the  mass  integral  in  F  igure  8  varied  hy  only  0.8'/  for 
the  second  through  fourth  stations,  and  the  momentum  integral  varied  hy  ♦  1  for  all  I  ,e  stations.  Note 
that  in  Figure  S  the  large  variation  in  the  integral  mass  balance  <  +  15.l'V )  at  the  initial  station  (/  =  1.4).  in 
part,  causes  both  £  and  <  to  diverge  from  the  consistent  results  obtained  at  the  other  stations.  The  degree  of 
consistent  exhibited  by  the  transport  coefficients  in  these  cases  is  particularly  impressive  when  it  is  noted 
that  the  scale  was  amplified  by  a  factor  ot  10  in  Figures  7  and  8  that  is.  the  magnitude  of  both  $  and  <■  is 
5  to  10  told  smaller  lor  these  cases  than  for  the  other  cases  in  Appendix  A.  in  which  jet  velocities  ranged 
from  050  to  over  3000  ft  sec.  (see  Table  l  >. 

A  survey  of  Figures  A- 1  through  A-10  shows  that  the  ranking  of  the  cases  and  the  total  percent  dev¬ 
iation  in  the  integral  mass  and  momentum  balances,  starting  with  those  exhibiting  the  most  consistent  (smooth! 
behavior  to  those  exhibiting  the  greatest  scatter  is  tabulated  here: 


SURVEY  OF  RESULTS 


Figure 

Maximum  %  Deviation  of 

1  (pUYAl. 
i 

Maximum  %  Deviation  of 

2  (P  U  U  Alj 
i 

A- 7 

12.5 

3.4 

A- 8 

4.0* 

• 

00 

rsi 

A  5 

16.4 

6.1 

A  2 

4.2 

3.3 

A  4 

29  7 

1.4 

A  10 

21.0 

4.0 

A- 3 

14.6 

7.1 

A  1 

28.3 

4.8 

A-9 

9.8 

1.4 

A3 

45.5 

7.7 

'Excluding  the  first  axial  station 


Scatter  in  the  last  four  entries  of  this  tabulation  is  great,  although  some  of  the  axial  variations  suggested  in 
the  corresponding  figures  may  be  real.  Generally,  the  less  deviation  in  the  integral  balances,  the  more  con¬ 
sistent  the  resulting  transport  coefficients.  Of  course,  the  axial  position  at  which  the  m;  imum  deviation 
occurs  is  also  important.  For  example,  a  large  deviation  at  an  intermediate  station  produ  *s  a  greater  scatter 
than  the  same  deviation  at  the  initial  or  final  stations  because  of  the  oscillatory  effect  on  the  axial  derivatives 
that  are  used  lor  caLulatum  of  the  Reynolds  transport  terms;  References  8  end  Additional  evaluation  of 
He  results  i<  required  be  tore  definite  conclusions  may  be  drawn:  however,  the  most  important  question  raised 
hy  these  results  is;  Will  the  simple  Composite  Trends  on  each  figure  permit  adequate  prediction  of  the  experi¬ 
mental  data  Ml  so.  the  detailed  variation  need  not  be  considered,  thereby  permitting  development  of  relatively 
simple  empirical  mixing  models. 

Ol  course,  the  validity  of  the  transport  coefficients  can  best  be  determined  by  using  them,  together  with 
initial  velocity,  concentration,  and  static  temperature  profiles,  to  predict  the  flow  field  from  which  they  were 
derived:  Reference  8.  Agreement  of  these  predicted  profiles  with  the  original  experimental  data  at  several 
axial  stations  is  conclusive  evidence  that  the  transport  coefficients  are  consistent  with  these  data.  Fstahlishim* 
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llu*  validity  <>|  a  vet  ot  transport  coefficients  in  this  manner  is  the  most  important  step  in  tire  development 
a'  iniXt 8  m<’do,s  Rosul,s  ol such  comparisons  lor  each  ol  the  cases  are  summarised  in  Tables  A-l  through 

xperimental  data  are  compared  with  the  predicted  <|uantities  at  various  axial  and  radial  positions  using 
the  Composite  Trends  of  Figure  1 3  lor  prediction  of  the  transport  coefficient*.  In  almost  all  cases  tire  initial 
'  .  U.  and  Fs  profiles  used  to  si;*  .  the  numerical  integration  were  the  first  axial  station  downstream  of  the 
core  region,  so  that  an  assessment  of  the  validity  of  the  £  and  <  initially  selected,  as  well  as  the  very  simple 
(omposite  Trends  ol  Figure  13.  could  he  marie  over  the  entire  transition  and  similarity  regions  (see  Figure 
-5Hor  which  experimental  data  were  available.  Results  presented  in  Tables  A-l  through  A-l 0  show  that 
the  (  omposite  Trends  were  indeed  adequate  for  the  prediction  of  both  Y  and  U  at  each  of  the  three  axial 
stations  tabulated,  which  included  the  last  statioi.  for  which  data  were  available.  Detailed  experimental 
temperature  measurements  were  not  obtained  so  temperature  comparisons  were  not  possible.  As  anticipated 
overall  agreement  between  experimental  an  I  predicted  Y’s  was  closer  than  between  experimental  and  pre¬ 
dicted  U  s  since  values  of  e  and  k/cp assumed  equal  to  £  were  used  in  the  calculations.  Nevertheless, 
adequate  overall  agreement  useful  for  hardware  predictions,  was  obtained.  Certainly.*  mixing  model 
applicable  to  the  range  experimental  conditions  presented  herein,  must  generally  predict  values  of  £  and  e 
of  the  magnitudes  of  the  Composite  Trends  of  Figure  1 3  (see  Section  li(-A-4|. 


The  C  omposite  £  s  plotted  on  the  “b”  figures  appear  to  be  an  adequate  representation  of  the  Com¬ 
posite  e  s  as  well 1  -  that  is,  an  Sc,  of  approximately  unity  is  the  best  estimate  that  can  be  obtained  from  these 
data,  except  in  Figure  5b  for  which  the  £  value  is  obviously  too  high.  Therefore.  Figures  5c  and  5d  were 
prepared  in  which  new  more  detailed  Composite  Trends  were  drawn,  using  the  same  points  as  Figures  5a 
and  5b.  Predictions  made  using  these  Trends  are  compared  with  the  experimental  data  in  Table  A-l  I .  Com- 
panson  ot  Tables  A-5  and  A-l  I  indicate  that  as  anticipated,  generally  better  agreement  was  attained  in 

h,  1  I',  °XLCP'  '°r  t"C1Y,s  prcdic,cd  at  the  z  =  5.37  in.  Predictions  at  this  station  were  somewhat  higher 

than  both  the  experimental  data  and  the  predictions  in  Table  A-5.  (Note,  the  mass  fraction  predictions  in 
a  e  ->  were  somewiat  low.)  Examination  ol  Figure  A-5c  shows  that  £  increases  modestly  with  axial 

*  u.  IOn  ,  ul,0.rlng  th,<i  c,fect  rc,;u,ted  ,n  ,he  Predictions  of  undermixing  in  Table  A-l  l ;  nevertheless,  results 
obtained  both  using  the  Composite  Trend  in  Figure  13  (Table  5)  and  the  more  detailed  Composite  Trends 
of  Figures  A-5c  and  A-5d  (Table  1 1 )  each  resulted  in  predictions  adequate  for  most  hardware  applications. 

mtv  .  Thc  re*u,,s  Panted  in  this  Appendix  conclusively  demonstrate  that  the  Inverse  Solution  Technique 

Z  k  f  SamPh"8  Pmbe  da,a  obtained  a«  a  few  a™«  nations,  as  long  as  the  integral  mass 

and  mom.-  ..urn  balances  at  each  station  are  reasonably  consistent.  However,  because  of  smoothing  used 
k  numerical  ditlerentiation  ot  the  data,  consistency  of  the  integral  balances  to  within  20'~,  appears 

adequate.  Ot  course  the  more  exact  these  balances,  the  more  details  of  the  behavior  of  the  transport 
coefficients  will  be  obtained. 
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Ikll  Aerospace  Company's  Advanced  Technology  Research  Denartmenf  hr.  ..  „„  , 

««'  ■  «" computations  cmOy:  te„«. ,  ^  T 

SiSSLr  ?“  ~  ■*  “•  — 1 «-  *~E22ZZ 

A.  COLD  FLOW  COMPUTATIONS 

Overall  Continuity 


dz 


d  id 

(pU>  +  -  —  l([,vn  =  0 

r  or 


Species  Continuity 


3Yi  dYi  I  a 

p  U  — +  pV  — L  =  1  1 
or  dr  r  dr 


3Yi  i  a  /  av,  i 

aT’7aT  —  j  » 


Axial  Momentum 
dU 


U  —  +  pV  — «  .  lit  +  i  _i  (  ,  .  .  dU  \ 

dz  3r  dz  r  dr  [  M  e)  ar  / 


P 

Energy 


(B-l } 


(B-2) 


( B-3 1 


♦rar  f  ,S,  wvf)  ♦  (w<*>p) 


(B-4* 
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with  miti.ii  and  boundary  conditions 


=  Utri.  Y, 

-  Yj  t r ►.  1  =  f  ir> 

(B-5;i| 

i.  ^  . 

dl 

dr 

dr  ' 

(B-5l<i 

=  Uc.  Y, 

=  «v,»c.  T  =  re 

cB  5ct 

Solutions  lor  the  velocity,  temperature  and.  mass  tractions  may  be  obtained  only 
alter  defining  the  molecular  and  turbulent  transport  coefficients  am*  *.ie  species  generation 
term.  u?r  file  molecular  transport  coefficients  of  mass*  momentum  and  thermal  energy, 
pl>j.  M.  and  kj.  respectively,  are  known  functions  of  composition  and  temperature:  whereas* 
their  respective  turbulent  counterparts  e  and  are  dependent  on  the  velocity  and 
density  fields,  initial  conditions,  and  system  geometry,  in  addition  to  the  composition  and  tem¬ 
perature  In  this  discussion  no  reactions  are  considered  anil  the  species  generation  term, 
c jj.  is  zero. 

One  of  the  dependent  variables*  V*  may  be  eliminated  from  the  above  equations  by 
transforming  to  the  von  Mises  coordinates*  (z,  i^),  defined  by 


,  ^  .* 

P  l  r 
dr 

V  —  =  P  v . 
dz 

d  d 

Substituting  for—  and  —  in  Eqs.  (B-2  to  B-4)usim* 
dz  dr 

d\  _  p  Vr  d\  +  d\ 

a//r.z  *  +  d7/0.i 


and 


\  P  Ur  4  \ 

3r/r./  V 


(B-(>a> 


(B-6h| 


<  B-“a ► 


<B-"h> 


II* 


defining  the  effective  viscosity*  /ieff  as 
^eff  =  M  +  f 


Bell  Aerospace  Company _ 

and  Hie  effective  turbulence  Prandtl.. Schmidt  and  Lewis  numbers 
O.  _  CP 


(k,  ♦  «,i 

"elf 

t&Vfrt 

Prj/Scj 

fiMter 

and  further  assuming  that  Sc,  and  Le,  are  independent  of  species  (therefore the 

Script,  i.  will  be  drop, led  for  these  terms  in  the  following  discussion). the soreraing  equations 
in  hie  i z\  p.<ine  are  given  by  ; 


Species  Continuity 

^i  j_  3_  T  Peff^Ur*  »Y,  *j  M.cjj 

3/-  &  b>fr  L  &  d\l/  J  +  p  u 

A*:ial  Momentum 

dU  =  I  dp  t  j_  3_  r  f*eff  p  U r*  au~| 
dr  pU  dz  ^  (_  *  d*J 

Energy 

c*  aT-  1  dP  |  1  3  [7  CP  \  ^eff  PUr:  QJ  1 

P  3z  p  dr  \l/  d\l>  L\Pr/  ^  a^J 

+  **eff p  Ur*  r/aun  +  _l_  3T  V  3Yi] 

**  |_W  /  s7  a 7  T  cPi 


fll-IOi 


tB-n  i 


t  Brill 


The  above  equations  ( B- 1 0  to  IM  2 .  were  soli ed  using  the  explicit  finite  difference 
senenu'  oi  Zieheie  and  Bleicl..  Ref  4K.  to  calculate  downstream  profiles  off.  Yj  and; 
r  given  initial  values  of  these  parameters.  In  this  formulation  the  difference  equations  fur 
each  dependent  variable.  F.  are  of  the  form 


Fj+ 1.  k  =  Am+  B(«»  Fj>k  +C<«>  FjJt+|  *  D<«>  Fj  k., 


(B-13 1 


I  iw 
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where  j  and  k  are  related  to  the  /  and  v  coordinates  respectively  tsee  sketch  below). 


U 


k  *  1 
k 

k  1 


1*1 


Sketch  .M  w.  /  Mesh  Network 

The  superscript.  S.  designates  the  specific  equation  considered,  i.e.,  «  -  1  for  species 

|IU 

continuity,  V  -  2  for  momentum  and  V  =  3  for  energy.  I  able  BI  lists  the  values  of  A 
and  l)^1  appropriate  for  each  equation,  where 

a  =  Mell  iJr3  * 

0  =  A  /  <k  (A^)3) 

On  the  axis  ot  symmetry,  v  *  0.  the  finite  difference  form  reduces  to 

,  U  >  „  *(V)  +  u  it)  t.  n+r  (C)  F  t 

!>.0  "  A0  *0  fj.0+co 

The  coefficients  Aq  are  identical  to  those  given  to  those  given  in  Table  B-l  whereas 
B0(t)  and  C  q  ( ^  are  given  to  Table  B-2. 

Transformation  from  the  ( z)  plane  to  the  physical  plane  is  obtained  using  the  quad¬ 
rature 


<  B-l  4a) 
(B-l  4b) 

(B-l  5) 


r  = 


r 

W  ‘  pu  J 


(B-l  6) 


The  governing  system  ot  equations  is  parabolic  hence  using  an  explicit  finite  dillvrcnci 
formulation  requires  that  a  stability  criterion  be  used  to  govern  allowable  grid  sizes.  A  /  and 
Aw.  Since  the  system  of  finite'  difference  equations  are  nonlinear  it  is  not  possible  to  deter¬ 
mine  by  analytical  means  the  precise  stability  requirements,  /leherg  and  Bleich.  Ret  48. 
obtained  an  estimate  of  the  criteria  tor  Az  and  Ai*/  based  on  linear  theory  and  numerical  experi¬ 
mentation.  They  found 

k  (A^)3 


Az< 


(-)  +  (-) 


(B-Pa) 
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and 


(A^)J  Sc 
Az<  - 7-  - 


(B-I7IU 


TA.H.I-:  H  I 

COEFFICIENTS  A‘*».  B<V).  C(V>  AND  D<V»  OF  EQUATION «B-I 

!  8=1  Species  Continuity  f 

A*  *  *  = 

St  Mj  cj,  i 

P  u  [ 

B(,)  » 

(S  jj  -(Sc  a)jt  k  +  a  -  (Sc  olj  k  .  J  ; 

j  C< ' »  * 

*i  (Sc  o)j  k  +  14 

|d*'»  « 

*#Sc  a)j.  k  -  '/i 

- — — - - 

8  =  2:  Axial  Momentum 

o 

II 

_  Az  /Jp\ 

'pU\dz/j+| 

b<~>  = 

C<2)  = 

S  °j.  k  +  Vi 

H 

>ia),k-  Vi 

8  '-3:  Energy 
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i  Mil  i  »•: 

< Oi  l  HCU  NTS  B(,<v  1  ANIH  „,v  ><>l  I  (R  ATION  <|M5) 


Equations  B  10.  B-I  I ,  and  B-!  2  are  solved  as  finite  difference  equations  (Lqs. 

(B-I 3>.  (B-I 5b  Tables  B-I  and  B-2)  numerically  on  the  digital  computer.  The 
program  requires  the  following  input:  Initial  velocity,  mass  fraction  and  total  temperature 
profiles,  static  pressure,  effective  Prandtl  and  Lewis  numbers  and  number  of  radial  mesh  points 
For  cases  where  1  e  =  0  small  but  finite  tree  stream  velocity  is  used  <  I of  maximum  velocity  L 
As  noted  above,  Zeiberg  and  Bleieh.  Reference  48.  obtained  an  estimate  for  the  step 
si/e  criteria  in  a\is>  mmetric  flows  based  on  linear  theory  and  numerical  experimentation.  How 
cut.  tins  step  si/e  criteria  required  that  Schmidt  number  j-.  be  less  than  or  equal  to  unity 
and  the  Lewis  number.  Le.  be  greater  than  or  equal  to  unity  \t  Bell  the  Step  si/e  criteria  was 
;jencr.iii/eO  to  appK  to  planar  and  a\ts>  nnuctne  geometries  !-t  arrMiary  values  of  Sc  and  be 
The  abilitx  to  voiKider  arbitrary  values  of  Sc  and  Le  is  important  in  studying  the  effects  of 
diffusion  coel t icienls  on  the  downstream  flow  field  de\el«  »p:nent .  and  was  required  in  the  de¬ 
termination  of  the  significance  ot  turbulent  transport  in  mixing  and  reacting  flows  presented 
in  Section  II I- A. 


S  A 
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B  REACTING  I  LOW  COMPUTATIONS 

I  lie  'lability  criterion  used  lo  determine  Hie  forward  marching  step  si/e.  U/j.  , .  does  not  con¬ 
sider  He  e  Meets  ol  reactions  in  its  calculation  Tins  . . .  occasionally  resulted  in  predictions  of  ,he 

mass  traction  being  negative,  winch  obviously  .s  not  physically  possible.  At  first,  this  difficulty  was  over¬ 
come  by  repeating  the  case  and  lightening  the  stability  criterion  in  an  arbitrary  manner.  This  practice 

frarn’niini Z  c°mpu,cr  "me  Therefore,  an  auxiliary  criterion  was  developed 

mum  "  I  ?  experimentation.  and  the  forward  step  si/e.  tAz  l.sc2.  was  expressed  as  a  function  of  mini- 
mum  velocity,  and  the  characteristic  ignition  time  ol  the  (low.  In  a  gaseous  reacting  How  with  typical 
scramjet  fuel  injection  conditions  the  ratio  (A/\, /t  has  been  found  to  be  as  high  as  13  0.  clearly 

illustrating  that  the  calculation  is  an  important  feature  in  the  solution  technique. 

The  Bell  computer  programs  arc  applicable  to  either  free  surface  or  wall  boundary  conditions  and 
a  transverse  coordinate  stretching  function  similar  to  (hat  used  by  Morgenthaler.  Reference  8.  has  been 
introduced  to  allow  liner  mesh  spacing*  in  regions  where  relatively  large  changes  in  the  dependent  variables 
occur  requiting  additionid  detail.  Note,  the  transverse  coordinate  stretching  capability  is  a  feature  not 
available  in  References  47  and  48.  7 

One  particular  aspect  of  the  Bell  effort  to  make  these  computer  codes  useful  engineering  design 

tools  is  the  procedure  tor  modeling  the  species  generation  term.  w.  when  finite  rate  kinetics  are  included 
A  limiting  ease  of  finite  rate  kinetics  occurs  when  the  reaction  rates  are  relatively  fast  so  that  thev  may  be 
assumed  to  be  infinite.  Ir.  this  case  equilibrium  conditions  may  be  assumed  which  require  a  simpler  and 
^expen^ecomputatton.  A  comparison  between  the  computer  cost  of  making  an  equilibrium  eomputa- 
'  "  "J  3.  ,n,tC  rjU  compu,at,on  (run  unt'l  equilibrium  conditions  were  reached )  shows  a  savings  factor 
of  80,.  |n  some  cases  even  an  equilibrium  computation  is  not  required  and  the  combustion  process  may 
be  idealized  by  considering  only  the  fuel,  oxidizer,  and  products.  Gosman.  et  aJ;  Reference  4©  Note 
however,  that  Spalding  recently  has  demonstrated  the  assumption  of  instantaneous  reaction  was  far  too 
crude  for  detailed  reacting  flow  analyses;  Reference  1 6. 

.  .  Wien  finite  rate  kinetics  are  used,  large  savings  in  computer  time  can  result  by  optimizing  the 
solution  for  the  species  generation  term,  since  it  must  be  obtained  at  each  mesh  point  in  a  reacting  tlow 
field  This  problem  may  be  restated  in  the  following  manner:  Find  the  most  efficient  solution  technique 
for  the  species  generation  term  given  the  chemical  reaction  steps  and  rate  constants.  Pergament  Reference 
51  has  shown  that  the  chemical  rate  equations  when  solved  explicitly  require  a  large  number  of  time  steps 
to  obtain  a  solution  due  to  the  “stiffness"  of  some  of  the  equations.  To  avoid  this  difficulty  the  implicit 
elution  technique  reported  by  Tyson  and  Kliegel.  Reference  52.  is  used  in  the  Bell  computer  code 
Bittkerand  Scullin.  Reference  53.  have  used  this  implicit  integration  scheme  in  developing  a  general 
chemical  kinetics  computer  program  where  only  the  reactions  and  rate  constants  need  to  be  specified  is 
input  to  define  a  'arge  number  of  gaseous  systems,  Efforts  are  in  progress  to  incorporate  this  generulitV 
into  the  Bell  mixing-reacting  computer  codes  so  lh.it  a  gre.it  curic-ty  of  reactions  may  be  considered. 

Hydrogen-oxygen  reacting  flows  are  of  particular  interest  in  a  number  of  air-breathim:  and  rocket 
propulsion  applications  that  is.  SCRAMJF.Tand  Attitude  Control  Systems.  Considerable  experience 
has  been  developed  at  Bell  in  optimizing  the  solution  technique  for  the  species  generation  term  of  this 
particular  gas  combination.  Table  B-3  lists  four  different  methods  used  to  solve  the  rate  equations  for  the 
H2/°2  reaction  sets  listed  in  Table  B-4A.  From  the  run  times  listed  in  Table  B-3.  it  is  clear  that  combining 
t he  q uasi-h neanzation  approach  suggested  by  Morretti.  Reference  54.  with  the  implicit  integration  method' 
leads  to  the  greatest  savings  in  run  time  -  for  example.  9.8  seconds  versus  the  604  seconds  required  when 
using  the  explicit  Runge-Kutta  integration  scheme  for  the  set  of  8  II-./0-.  reactions  considered 
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TABLI  B-3 

COMPARISON  Ol  KIM  I  K  COMIHTATIONAL  1 1  ( IINIOt'l  S 


Description  of  Solution  Technique 

Computer  Run 
Time  (sec) 

1.  Quasilineart/ed  Rate  Equations  w*th  Implicit 

Integration  Scheme 

9  8 

2.  Quasilinean/ed  Rdte  Equations  w«th  Exact  Solution  | 

horn  E igenvatue Solution 

45.1* 

3.  Complete  Rate  Equations  (nonlinear)  With  Implicit 
Integration  Scheme 

1190* 

4.  Complete  Rate  Equations  (nonlinear)  with  Explicit 
Integration  Scheme 

604.0 

^Computer  programs  using  these  solution  techniques  perform  auxiliary 
computations  that  are  not  required.  Removal  of  these  auxiliary  com¬ 
putations  would  in  the  trest  case  reduce  computer  run  time  by  50%. 


TABLF  B-4A 
RATH  CONSTANTS 
IK/AIR  BAC  SURVHY 


Reaction 

Number 

Reaction  Rate  Parameter* 

<CGS  Units) 

Reaction 

A 

N 

E 

1 

H* 

■f 

°2 

-  OH 

+ 

° 

2.135-1014 

0.0 

16,533 

2 

0 

*■ 

H2 

*  OH 

+ 

H 

2.512-1012 

0.0 

7,672 

3 

OH 

+ 

h2 

*  h2o 

H 

6.310-10’3 

0  0 

5,870 

4 

OH 

+ 

OH 

=  h2o 

+ 

0 

7.575-1012 

00 

996 

5 

M 

+ 

H2 

*  H 

+■ 

H 

2.202-1019 

•1.00 

106,128 

6 

M 

+ 

h2o 

=  H 

+ 

OH 

7.690-1020 

1.00 

121.770 

7 

M 

+• 

OH 

=  H 

+■ 

0 

6.700-1018 

*1.00 

104,742 

8 

M 

°2 

=  0 

0 

1.000-1020 

•1.00 

121,374 

Forward  Rate 


ATN  exp  (-E/RT) 


Reverse  Rate 


Forward  Rate 


Equilibrium  Constant 


A  study  also  was  made  to  optimize  the  iteration  method  for  the  density  at  the  succeeding  step, 
which  is  required  as  a  consequence  of  the  quasi-lineari/ation.  It  was  found  that  usually  a  maximum  of 
two  iterations  were  required  for  convergence  (temperature  tolerance  of  IK)  when  using  the  method  of 
false  position.  However,  frequently  no  iterations  at  all  were  required  when  the  value  of  the  previous  step 
was  used  as  the  first  estimate. 

Of  course  from  a  physical  standpoint,  one  of  the  most  difficult  problems  in  defining  the  species  gen¬ 
eration  term  is  specifying  the  detailed  chemical  kinetics  that  is.  the  reaction  steps  and  their  rate  con¬ 
stants.  Since  the  kinetics  determine  the  history  of  the  gas  composition  and  temperature,  a  study  was 
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made  to  determine  the  best  set  ol*  rate  constants  consistent  with  the  eight  reactions  used  in  the  quasi- 
InuMri/cil  approach.  As  a  basis,  the  I H  reactions  and  corresponding  rate  constants  reported  by  Baulch. 
el  al.  Reference  55  through  58.  Table  B4B,  and  recommended  by  Biftker.  Reference  53.  as  the  most 
appropriate  set  of  If,  '<),  reactions,  were  used  to  obtain  the  time  history  of  the  reaction.  Predictions 
wen'  also  made  using  the  eight  reactions  and  rate  constants  which  have  been  reported  in  the  analysis  of 
hydrogen-air  Hows.  A  set  of  rate  constants  designated  as  Bell  Survey.  Table  B4A.  which  gave  reasonable 
agreement  with  the  predictions  obtained  using  the  18  reactions  (see  Figure  B-l )  were  thus  determined 
Further  comparison  showed  that  the  set  of  rate  constants  proposed  by  Fern,  et  al.  Reference  59.  Table 
B4C  predicted  the  time  tor  the  reaction  to  reach  equilibrium  was  270  microseconds  sooner  than  the 
values  obtained  using  the  1 8  reactions  (see  Figure  B-l  >.  Note,  that  t his  difference  corresponds  to  a  dis¬ 
tance  of  0.54  ft  in  a  SCRAMJFT  combustor  with  a  velocity  of  2000  ft /sec. 


TABLF  B4B 
RATE  CONSTANTS 

H2/A1R  -  BITTKER  AND  SCULLIN,  REF.  S3 


Reaction 

Number 

Reaction 

RMCtKHi  Raw  Paramawrt 
(CGS  Unittl 

A 

N 

E 

t 

Hz 

OH 

*  h20 

4- 

H 

2.10-10>3 

f BM 

mgm 

2 

H 

°2 

-  OH 

4- 

0 

1.25*1014 

3 

0 

H2 

OH 

4- 

H 

2.96-1013 

Bill 

9,800 

4 

H 

+ 

°2 

-  ho2 

4- 

M 

1.59-1015 

-1.000 

5 

H 

♦ 

H 

11 

I 

N) 

4- 

M 

6.52-1023 

3.770 

6 

H2 

+ 

ho2 

*  H2°2 

4* 

H 

960-1012 

24,000 

7 

M 

+ 

=  OH 

4- 

OH 

1.17-101^ 

0.0 

45.500 

8 

H 

4- 

ho2 

-  OH 

4- 

OH 

7.oo-io!3 

0.0 

0 

9 

H 

4* 

OH 

-  h2o 

4 

M 

7  50-1023 

2.60 

0 

10 

0 

+ 

0 

-  °2 

4* 

M 

1.38-1018 

1.00 

340 

11 

0 

+ 

H2° 

-  OH 

4* 

OH 

5.75-1013 

0.0 

18,000 

12 

H2 

4* 

°2 

=  OH 

4* 

OH 

1.00-1013 

0.0 

43.000 

13 

OH 

+ 

ho2 

*  H20 

4* 

°2 

6.00-1012 

0.0 

0 

14 

0 

+ 

CM 

o 

X 

*  OH 

4 

°2 

6.00-1012 

0.0 

0 

15 

ho2 

♦ 

ho2 

=  H2°2 

4- 

°2 

1.80-1012 

0.0 

0 

16 

OH 

4- 

H2°2 

-  H20 

4* 

ho2 

1  00-10]3 

0.0 

1.800 

17 

0 

+ 

H2°2 

=  OH 

4- 

ho2 

8.00- 10>3 

0.0 

1.000 

18 

H 

+ 

HjOj 

-  h2o 

4- 

OH 

3.18-IQl4 

0.0 

9.000 

Forward  Rate  =  ATN  exp  (-£  RT) 


Reverse  Rate 


Forward  Rate 
Equilibrium  Constant 
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TABLE  B-4C 
KA  I  F  CONSTANTS 
Hv AIR  •  W.KKI.etal  RLF59 


Reaction 

Number 

Reaction 

Reaction  Rat*  Parameter! 

(CCS  Units) 

A 

N 

E 

1 

H 

f 

°2 

=  OH 

♦ 

0 

3  000*1014 

0.0 

1  7,443 

2 

0 

f 

H2 

-  OH 

4 

H 

3. 000* 10 >4 

0.0 

7,979  I 

1  3 

OH 

f 

H2 

-  H20 

F 

H 

3000*I0<4 

0.0 

5,979  * 

4 

OH 

+ 

OH 

=  H20 

F 

0 

3.000M014 

0.0 

5,979 

5 

M 

f 

H2 

*  H 

F 

H 

1. 850*1020 

too 

106,920 

6 

M 

4- 

h2 0 

»  H 

■f 

OH 

9.660M021 

-too 

123,156 

7 

M 

♦ 

OH 

*  H 

■F 

0 

8  000M019 

-t.oa 

(03,356 

8 

M 

♦ 

°2 

*  0 

F 

0 

5.800*  I0>9 

T.OO 

(19.988 

Forward  Rate 


Reverse  Rate 


ATN  exp(-E/RT) 

Forward  Rate 
Equilibrium  Constant 


,  ,  made  for  tw>th  mixing-reacting  and  frozen  hydrogen/air  flows  using  the  computer 

code  described  above  are  shown  in  Figure  B-2  to  illustrate  this  computation  capability.  The  initial  hydro 

dlme^in  mi  rTf  U"d  T  tha'  rcp°r,ed  by  Ropcrs-  R*fere"«  60,  on  the  center  plane  of  a  three- 
dimensional  flow  field  at  seven  diameters  downstream  of  the  injection  of  hydrogen  into  a  supersonic  air 

stream.  Assuming  this  same  profile  was  the  initial  profile  for  a  two-dimensional  fiat  plate  with  a  trans¬ 
verse  injection  slot,  predictions  of  the  hydrogen  mass  fraction  at  30  diameters  downstream  were  made 
assuming  both  Irozen  or  reacting  How  conditions.  Ignition  and  combustion  were  achieved  by  specifying 
a  uniform  total  temperature  of  2000°K  across  the  boundary  layer.  Of  course,  in  actual  rocket  engines 
ignition  is  accomplished,  in  part,  by  entrainment  of  the  hot  combustion  products  in  local  recirculation 
regions.  Consideration  of  this  effect  is  the  subject  of  a  current  Bell  investigation. 

Finally,  it  should  be  noted  that  the  Bell  computer  codes  contain  a  number  of  useful  auxiliary  sub¬ 
routines  These  routines  provide  useful  engineering  information  although  they  do  not  play  a  role  in  nre- 
dieting  the  now  field  development.  Two  specific  examples  are: 


•  The  calculation  of  combustion  effeciency  shown  in  Figure  B-3  using  the  Combustor  Cor¬ 
relation  Technique,  by  Morgenthaler.  et  al.  Reference  19.  which  assumes  cold  nonreactive 
mixing  of  fuel  and  oxidizer  to  be  the  controlling  process  in  combustors  empfovine  diffusion 
flames  and  *  - 


•  The  calculation  of  entrainment  rates  in  heated  and  cold  jets  shown  in  Figure  B4  reported 
by  Pescli ke :  Reference  61, 
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Hydrogen  Mass  Fraction 
<%> 


f  igure  B-2.  Predicted  Hydrogen  Mass  Fraction  Distributions  at  /  D  =  30,  Using  Experimental 
Center  Plane  Initial  Data  of  Rogers,  Ref  60. 
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M’l'IADIV  ( 

Sl1 1*1*0  R  flN(i  l'K(X.KAM.  I  M'lKIMI  MAI.  SIT.'DV  OF  PARTICLF  DYNAMIC'S) 


In  .ni  experimental  program  Bell  is  study mg  Moutange,ii,.i|  micclion  «t  tw.vpl.ase  and  single  pl.a 
K..  .nr.;  subsonic  llou  |  h.s  mveslig.it  ion  .inns  to  tin.l  mlormalion  lor  the  rational  design  of  combustion 

vlKlHHVts  dtul  HI  KM  (Ms  U\IM|£  powdered  llit/iv 


I  igiiret-l  shows  photograph  ol  the  experimental  tacil.lv  which  has  just  been  completed  A  centn- 
lugal  blower  puts  air  in  a  4.25-m.  square  test  section  with  glass  windows  on  all  four  sides.  So  it  s  possible 
to  make  simultaneous  observations  at  right  angles  to  each  other  by  regular  or  schheren  photography  Figure 
(  shows  the  double-view  schhetvn  system.  It  has  the  elements  of  a  standard  schheren  system  but  two 
plane  mirrors  have  been  added.  I- me  adjustment  of  the  plane  mirror  No.  2  makes  it  possible  to  use  the 
same  km  c  edge  lor  both  views;  Reference  (-2.  Removal  of  the  knife  edge  converts  tbissvstem  to  one  for 
regular  photography  wild  two  views. 

Traverses  across  the  test  section  with  a  pitot-static  tube  indicated  satisfactory  velocity  distribution 
across  the  test  section  as  shown  by  Figure  (  '-3.  Reference  values  indicated  in  the  figure  were: obtained  from 
a  stationary  upstream  pitot-static  probe  at  the  venter  ot  tic  cross  section.  Deviations  of  *  2'>'  from  the  refer¬ 
ence  value  are  also  entered  in  the  figure  Tins  equipment  is  being  used  to  observe  the  jet  and  its  spreading  m 
two  dimensions.  I  he  jet  behavior  will  be  studied  for  various  parameters  that  influence  the  trajectory,  such 
as  stream  and  jet  velocities,  particle  loading,  size  and  concentration,  as  well  as  the  injection  angle  Figure 
-4  shows  a  preliminary  schheren  photograph  of  a  mitogen  let  carrying  glass  beads  I -325  mesh  i  at  a  mass  ratio 
ot  about  |()0;  I .  the  jet  velocity  ,s  about  HO  ft  s  and  the  air  velocity  VO  It  s.  Views  at  right  amrles  ,  Figure  C  -2 1 
will  he  obtained  alter  installation  of  high-quality  plane  mirrors 

One  ot  tiie  lirst  experiments  w  ill  be  the  observation  of  two-dimensional  jets.  Figure  C-5  shows 

schheren  photographs  ol  a  helium  jet  o|  I4~(J  ft  s  into  different  air  Hows  (36  to  250  ft  si  where  the  w  idth 
ot  the  injection  slit  is  0  0095  in. 

Powdered  fuels  must  he  efficiently  packed  into  a  fuel  tank  We  made  attempts  to  pack  boron 
powders.  Results  with  amorphous  boron  indicate  that  this  material  forms  large,  low-density  agglomerates 
which  lead  to  packing-density  ratios  (weight  of  packed  sample  weight  of  same  volume  of  solid  material  i  of 
less  than  40  ,  even  with  vibration  or  surface  treatment.  C  rystalline  powder  of  100  mesh  gave  a  density  ratio 
ol  a  tout  3  -  ter  surface  treatment  while  a  325  mesh  sample  yielded  only  56'V.  Mixtures  of  these  two 

samples  did  not  indicate  a  composition  of  more  efficient  packing  readily  obtainable  with  other  materials.  A 
report  on  this  program  has  been  prepared:  Reference 

Spheroidi/ing  of  particles  improves  the  packing  efficiency  of  other  materials  and  might  also  be  bene¬ 
ficial  or  boron  However,  boron  and  bo, on  compounds  are  extreme, y  ditf.cult  to  spheroid, /e  becanv  ol 
their  hicn  melting  point  and  low  .tens, tv  I  tie  o,  producing  small  samples  for  evaluate.,  was  pr..i„b,me 
lot  tile  present  program.  Finally,  alter  completion  ot  this  program,  two  small  samples  of  spheroidi/ed  boron 
were  received  that  had  been  produced  previously  by  TAI  \  Division.  1 1  urn  phrevs  Corporation.  In.e  maximum 
partite  si/e  was  about  s  urn  lor  both  sample-,  with  average  si/es  of  about  20  and  40  jim.  Packing  ol  these 
sampies  yielded  a  packing  efficiency  of  (,~  ■  for  the  larger  particles  and  64  for  the  smaller  ones  \-  j„  ,|K. 
case  ol  crystalline  boron,  the  larger  average  particle  size  y , elded  better  packing  efficiency.  Although  these 
packing  etlic.enc.es  are  no  better  than  those  for  crystalline  boron,  spheroid, zing  may  well' have  improved  the 
packing  efficiency  beyond  that  obtainable  with  the  original  material.  Unfortunately,  no  information  con¬ 
cerning  the  original  boron  could  be  obtained  except  that  ,t  had  been  produced  some  time  ago  as  P  .rt  of  , 
urge -nc ale  commcrcidl  production 
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DOUBLE  -  VIEW  SCHUEREN  SYSTEM 


Figure  C-*.  Double-View  Schlieren  System 
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Distance  from  Centerline,  in 

Figure  f-3.  Velocity  Distribution  in  t !>e  Test  Section 
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R«"*C4.  Sc*l**mn  of  «  Jet  of  (3—  Beadi Canted  by  hfitrogen (Mm  Ratio 

^Wwtaitily  100:1) 
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(b)  150  ft/s 


CO  250  ft/s 


FigU»C-5.  Schlieren  Photographs  of  a  Helium  Jet  ( 1470  ft/s)  into  Air  Flows 
(the  width  of  the  injection  slit  is  0.0005  in.) 
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